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PROGRAM  IA  -  COMPUTATION  OF  DOWNWASH  DUE  TO  TRAILING  WAKE 


This  program  computes  the  downwash  at  the  rotor  blade  using  Equation  15  of 
Ref.  1  ?  In  this  program  the  symbol,  Vj,  is  used  for  Wj  of  Ref.  1  and  x  for  ^  . 

The  chord  variation,  x,  of  Ref.  1  is  here  assumed  to  be  zero.  The  computational 
techniques  involved  in  the  solution  of  this  equation  by  numerical  methods  are  chiefly 
those  of  integration  and  harmonic  analysis.  It  was  decided  to  perform  the  numerical 
integration  by  Simpson's  rule  (an  approximation  of  the  curve  by  a  series  of  second  degree 
parabolas)  since  this  method  seemed  most  appropriate  to  the  nature  of  the  function.  Ex¬ 
pressed  symbolically,  this  means  that  given  a  function,  y  =  f  (x),  then 


r 


y  d»  =  y  [yD  +  yn 


n  -  1 

m2 

i=  1 


n  -2 

•>£  v 


i  =  2 


where  yQ  =  f(a)  and  y^  =  f(b),  i  =  odd  integer,  j  =  even  integer 

For  ease  in  programming,  a  separate  subroutine  was  written  to  perform  the  harmonic 
analysis.  This  was  done  as  follows: 


Let  g(x)  represent  the  function  to  be  analyzed, 

then  g(x)  =  a^  (a^  cos  k  x  +  bj^  sin  k  x)  where  the  a^ 

and  b^  are  defined  in  the  following  relations: 


Q  +  2ir 


°o  2h 


g(x)  dx 


Q  +  2ir 


g(x)  cos  k  x  dx 


g(x)  sin  k  x  dx 


(k/0) 


£  R.  H.  Miller.  "Theoretical  Determination  of  Rotor  Blade  Harmonic  Airloads.  " 
M.  I.T.  ASRLTR  107-2,  1964. 
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This  analysis  Is,  of  course,  applicable  only  for  a  range  of  2w  (Q  Q  +  2tt).  The  Inte¬ 
gration  again  was  performed  using  Simpson's  rule. 

The  programs  are  written  in  Fortran  II  and  should  be  operable  at  any  IBM  709  or 
7090  installation  with  one  exception.  The  Massachusetts  Institute  of  Technology  Compu¬ 
tation  Center  at  which  these  programs  were  developed  assigns  logical  tape  numbers  to 
physical  tapes  as  follows: 


Physical 


Logical 


1 

4 
2 
8 

5 

6 
7 
3 


input 

output  for  printing 


output  for  punching 


Interval  sizes  of  7.  5°  for  ^  and  7. 5°  and  2.  5°  for 0  used  with  an  m  of  3  and 
QT  of  -  have  produced  results  of  sufficient  accuracy  for  the  present  purposes. 


CAUTION:  Care  must  be  taken  in  selecting  interval  sizes  so  that  an  odd  number  o*  points 
is  to  be  used  in  each  integration. 


The  Function  V,  (y^x)  ~  Trailing  Wak< 


r>2  tt  m  +  y  [  t  +  d  cos  <P  ~  x  cos  (Jc  -  (p)  +  p  sin  ]  -  px  smif}  $ 

^  f  .  2  ,2  2  _  0  Q  x  imO. —  _ ia 2 


[l  +x  +d  +  z  ~2x  Q  +  2d(icos^>-x  cos^r)] 
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where  Kj  =  Kc  cos  n (j>  +  K$  sin  n <j) 
d  =  p  [  (2it  m  +  yj/'  +  y)  -  <p] 

z  =  Xi  (2wm  +  ^  +  y)  -  ^  ]  (1 . 0  +  a  ccs^ )  -  bQ  (/-  x  ) 

This  program  computes  the  value  of  the  function  for  Kc=  r-  I  for  each^  -  x  pair  used  at  from 
one  to  four  values  of  y  for  a  given  number  of  values  of  j/C  A  value  of  zero  is  assumed 
to  be  among  the  y's  read  into  the  computer  for  each  £~x  pair,  so  that  if  only  one  value 
of  y  is  to  be  used  that  y  must  be  zero.  The  integration  for  the  y  =  0  case  only  is  split 
into  two  sections  and  integration  performed  between  limits  of  j/f+y— >2trm+yr+-y  +  QT 
and  2nm  +  j/r  +  y  +  QT  — *  2ttm  +  +  y.  The  results  of  the  latter  integration  are 

termed  ZAPSI  (zq(^))  for  the  cos  component  and  ZBPSI  for  the  sin  component. 

The  results  of  integration  for  the  first  section  and  the  results  of  integration  for  any  other 
y's  are  termed  YiPSI  (i  =  1 , 2, 3,4, 5,6, 7, 8).  The  odd  integers  are  assigned  to  the  cos 
components  and  the  even  integers  to  the  sin  components.  The  i's  will  correspond  to  the 
y's  in  the  order  in  which  the  y's  were  read  into  the  computer,  i.e. ,  YIPSI  and  Y2PSI  will 
be  assigned  to  the  results  for  the  first  y  read  into  the  computer. 


Next  the  YiPSI  terms  are  summed  at  each^f  .  The  sum  for  the  cos  components  is 
called  CHI  )  and  the  sum  for  the  sin  components  is  called  SIGMA  (o’).  Harmonic 
analysis  is  then  performed  on  each  of  the  four  functions,  ^  ,  cr,  and 

The  time  required  to  obtain  results  from  the  IBM  7090  computer  for  one  Jt — x 
combination  with  3  y's,  m  =  3  and  interval  sizes  of  7.  5°,  7.  5°  and  2.  5°  is  approximately 
two  minutes. 


Subroutines  required  with  this  program  will  be  described  in  the  following  pages. 

The  listing  of  the  program  in  Fortran  II  >gether  with  the  subroutines  is  given  in 
Appendix  A. 

Typical  results  are  listed  in  Appendix  B. 
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Phyf  1  cal  Explanation  of  Symbols 


Y  &  rotor  azimuth  measured  from  blade  downwind  position 

y  blade  spacing  (cf.  J  ) 

m  number  of  wake  spirals 

V,(^x  )»w,  (^,») 
jO  rotor  span  parameter 

^  azimuth  of  wake  measured  from  downwind  position 
x  23^  rotor  span  parameter  (y  of  Ref.  1) 

n  25J  harmonic  of  rotor  speed 

p  252  advance  ratio 


a  /v< 


z 


J 


b 

o 


*0 

vertical  distance  travelled  by  rota  hub 
horizontal  distance  travelled  by  rotor  hub 
coning  angle 
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Input  Format 


Card  No. 

Card  No. 
Card  No. 

Card  No. 

Card  No. 

Card  No. 

Card  No. 

Card  No. 


1  L  KL  NP  II 

12  12  12  12 

1A  BO  (FORMAT  E8.  3) 

2  EMU  AMBDA  EM  DEPSI 

E8.3  E  8. 3  E  8. 3  E13.8 


MM 

12 

DELPHI 
E  13.  8 


NN  QT 

12  E  13.  8 

DELTA  K  IJ 
E 13. 8  13  12 


3  A  IK 
E8.3  12 


4  ELK  IL 

E8.3  12 

5  EN  IM 

E8.3  12 

6  X  IN 

E8.3  12 


7  GAMMA 
E  13.  8 
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Explonation  of  Symbols 


L 

KL 

NP 


MM 

NN 

QT 

BO 

EMU 

AMBDA 

EM 

DEPSI 

DELPHI 

DELTA 

K 

IJ 


number  of  harmonics  to  be  computed  in  the  Fourier  analysis 

+  1  — >  print  optional  table  No.  2,  -1  — >  suppress  this  print  out 

interval  at  which  table  No.  2  is  printed  (i.  e.  ,  +  1  will  cause  printout  of 
integral  at  every  Iff  used,  +  2  at  every  other  ~jlr  etc.) 

number  of  times  card  No  2  is  to  be  repeated 

+  1  — ^  print  optional  table  No.  1,  -  1  — *  suppress  this  print  out 

interval  at  which  table  No.  1  is  printed  (see  NP) 

determines  limits  of  integration  for  y  =  0  case  (i.  e.  ,  QT  =  -  ^  will  cause 
integration  from  )lr  +  y  to  y  +  2irm  -  ir/2  and  from  7 lf+  y  +  2irm  -  rr/2 
to  y/  +  y  +  2  tt  m) 

b 

o 

M 

X 

m 

interval  size  used  for  ^ 

interval  size  used  for  p  [  in  y  =  0  case,  in  region  ^  +  y  -> ^  +  y  +  2 it  m  +  QT] 

interval  size  used  for  ©  for  y  =  0  case  in  region  0r+y  +  2iTm  +  QT  to 
+  y  +  2  it  m 

number  of  points  of  ijf  to  be  used  in  harmonic  analysis 
number  of  values  of  a  to  follow 


A 

IK 


number  of  values  of  £  to  follow 


ELK 

IL 


i 

number  of  values  of  n  to  follow 


EN 

IM 


number  of  values  of  x  to  follow 


X 

IN 


number  of  values  of  y  to  follow 


GAMMA 
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Output  Formot 


MU  =  LAMBDA  =  M  = 

DELTA  PSI  =  DELTA  PHI  (1)  = 

CHI 

A  -  ZERO  = 

A  - k=  B  -  k  = 

I  1 


SIGMA 
A  -  ZERO  = 
A  -  k  = 


1 


ZAPSI 

A  -  ZERO  = 
A  -  k  = 


i 


i 


ZBPSI 

A  -  ZERO  = 
A  -  k  = 


I 


i 


(Optional  Table 


CHI  SIGMA 

i  L 


(Optional  Table  No.  ^ 
PSI  Y1PSI 


1 


I 


PSI 

i 


Y2PSI 


N  =  L  =  A  =  X 

DELTA  PHI  (2)  = 


ZAPSI 


1 

Y3PSI 

i 


ZBPSI 

I 


Y4PSI  etc 

1 
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DELTA  PHI  (1)  »  DELPHI 
DELTA  PHI  (2)  ^  DELTA 


[  see  explanation  of  symbols  (input)  ) 


A  -  ZERO  ;t#  constant  term  resulting  from  Fourier  analysis 
A  -  k  coefficients  of  the  cos  k  terms  resulting  from  the  Fourier  analysis 

B  ~  k  «  coefficients  of  the  sin  k  "\lr  terms  resulting  from  the  Fourier  analysis 
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Subroutine  Harnal  (PSI,  FPSI) 


This  subroutine  performs  the  actual  harmonic  analysis  as  previously  described. 


Information  is  transmitted  to  and  from  this  subroutine  by  means  of  a  common  state¬ 


ment. 


where 


COMMON  PSI,  FPSI,  I,  L,  K,  DEPSI,  AQUAY,  BQUAY,  A-ZERO 

PSI  ~  variable 

FPSI  -  a  function  of  PSI 

I  -  an  indicator  to  shc*v  type  of  function  so  that  excess  computation  is 
avoided,  i.e.,  +1  — *  an  even  function,  ♦  2  — ►  an  odd  function, 

♦  3  function  is  neither  even  nor  odd.  This  is  automatically  set  in 
(MAIN) 

L  _  number  of  harmonics  to  be  computed  '  (maximum  =  20) 

K  ~  number  of  values  of  PSI  to  be  used,  must  be  an  odd  integer  >  3 
DEPSI  -  interval  size  used  for  ~jfs 

AQUAY  -  output  -  the  coefficient*  a-k  of  the  cos  k  terms 
BQUAY  -  output  -  the  coefficients  b~k  of  the  sin  k^r  terms 
A-ZERO  -  output  -  the  constant  term 


It  is  assumed  that  the  table  of  wi| describe  a  range  of  2w  if  the  function  is  neither  even 
nor  odd  and  otherwise  a  range  of  w. 
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Subroutine  Prince  (L.  AQUAY.  BQUAY.  A  ZERO 


This  subroutine  prints  out  the  results  of  the  harmonic  analysis  in  the  following  form 


A-ZERO  = 


A  -  1  = 


constant  term) 


B  -  1  = 


A  -  2  = 


B  -  2  = 


A  -  L  = 


B  -  L  = 


(coefficients  of  cos  k  ~jf  terms)  (coefficients  of  sin  k  ~jf  terms) 


An  optional  version  also  punches  these  results. 


Subroutine  TRIP  (CHI#  SIGMA,  PSI,  ZAPSI,  ZBPSI,  NN,  K) 


This  subroutine  prints  optional  table  No.  1. 


Subroutine  PRICK  (NP) 

Information  is  transmitted  to  this  subroutine  via  COMMON.  It  is  used  to  print  op¬ 
tional  table  No.  2.  The  subroutine  must  conform  to  the  number  of  y's  used  in  MAIN.  If 
3  y's  are  used  PRICK  must  output  PSI  and  ViPSI  (i  =  1  to  6).  The  present  version  is  written 
fcr  three  y's,  hence  some  modification  will  be  necessary  if  a  different  number  of  y's  is  to 
be  used. 
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Somple  Data  for  Program  IA 


Card  No. 
Card  No- 
Card  No. 

Card  No. 
Card  No. 
Card  No. 
Card  No. 
Card  No. 
Card  No. 
Card  No. 
Card  No. 

Card  No. 
Card  No. 
Card  No. 
Card  No. 
Card  No. 
Card  No. 
Card  No. 

Note  tfiat 
(cards  No. 


1  10+1  +  1  +  2  +  1  +  1-  .  15707963  E  +  1 
Ifl  -f .  Ofc9  E  +  o 

2  +  .  100  E  +  0  +  .  500  E  -  1  +  .  300  E  +  1  +  .  13089969  E  +  0  +  .  13089969  E 
+  0  +  .  43633233  E  -1+49+1 

3  +  .  500  6  +  0  +  1 

4  +  .  100  E  +  1  +  1 

5  +  .  300  E  +  1  +  1 

6  +  .  550  E  +  0  +  3 

7  +  .  0000  0000  E  +  0 

8  +  .  20943950  E  +  1 

9  +  .41887901  E  +  1 

10  +  .  100  E  +  0  +.200  E  -  1  +  .  300  E  +  1  +  .  13089969  E  +  0  +  .  13089969  E 
+  0  +  .43633233  E  -  1  +  49  +  1 

11  +  .  000  E  +  0  +  1 

12  +  .  100  E  +0+  1 

13  +  .000  E  +  0  +  1 

14  +.975  E  +  0  +  3 

15  +.  0000  0000  E  +  0 

16  +  .  20943950  E  +  1 

17  +  .41887901  E  +  1 

the  data  block  (cards  No.  3"+  9  )  following  card  No.  2  must  be  repeated 
11— *17  )  even  though  some  of  the  parameters  are  unchanged. 
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PROGRAM  IB  -  COMPUTATION  OF  DOWNWASH  DUE  TO  SHED  WAKE 


This  program  computes  the  downwash  at  the  rotor  blade  using  Equation  21  of 
Ref.  1  with  the  symbol,  V substituted  for  and  x  for  as  before. 


y  +  2  it  m 

V2  "  /  K2 KVm(^ \\ “  "7 " A'T 

'  I  z  *-  (x  sin \y-<P)  *  d  siny>) 


f 


) 


f~  2  +  d  cos£ 

^  -  X  COS  {ft~ft) 

d  cos  ft  ~  x  cos  (V'"-^)] 

+d  -2x  ^  coi{ft-ft)  +  2d(^  cos^>  -  x  cosi 

h  | 

/  2  2  ,2  0  ,  , 
x  *-z  •‘•d  ~2xd  cos TJ 

!r 

d  ^ 


d  =  p  [(2  tt  m  ^  +  y)  -  (j)  ] 

z  =  ^[(2rrm^yr  +  y)  -ft] 

=  cos  n^>  +  sin  n ft  (K  and  K  are  set  to  1.0) 


This  program  evaluates  the  integral  in  two  steps.  For  the  first  step  (from  f  y 
2TTm+yf+y  +  QT  )  integration  Is  performed  by  Simpson's  rule  using  <2  a*  defined  above. 
The  results  of  this  integration  are  labelled  S 1  PS  I  for  the  cos  component  and  S2PSI  for  the 
sin  component.  For  the  second  step  cos  n  ft  and  sin  n  ft)  are  expanded  to 

sin  nlTp  -ft)  sin  n  ft  +•  cos  n  {ft -ft)  cos  n ft  and 

cos  n  (iff  -ft )  sin  n  ft  -  sin  n  {ft -ft)  cos  r.lft  respectively. 

The  integration  is  then  carried  out  from  (ijr  +  y  +  2  v  m  +  OT)  to  ( ft  t  y  t  2  u  m)  using 

=  Kc  sin  n  {ft  -ft)  siring  +  (  -K  J  sinn  {ft -ft)  cos  n  ft .  The  coefficient  of  Kc  is  cal  led 
S4PSI  and  the  coefficient  of  ("Kj)  is  called  S3PSI.  The  value  of  the  function  at  the  upper 
limit  is  approximated  by  setting  ft)  =  ft  -  1  2  A  ft).  This  gives  a  very  good  approximation 
when  the  interval  size  used  in  this  region  is  on  the  order  of  2.  5°.  Harmonic  analysis  is  then 
performed  on  all  fo^r  of  these  functions  (S1PSI,  S2PSI,  S3PSI,  S4PSI).  This  program  does 
not  consider  those  terms  involving  cos  n  ( ft  -  ft). 

The  time  required  for  one  case  on  the  IBM  7090  computer  (one  x,  one  2,  one  y)  is 
approximately  .  5  minutes  for  a  7.  5°  interval  size  and  m  =  3. 
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Subroutines  HARNAland  PRINCE  are  used  with  this  program.  The  program  used 
for  Vj  is  listed  in  Appendix  C  with  typical  results  in  Appendix  D. 


Input  Format 


Card  No.  1 

L 

NP 

KL 

II 

QT 

12 

12 

12 

12 

E  13.8 

Card  No.  2 

EMU 

AMBDA 

EM 

DEPSI 

DELPHI 

DELTA 

K 

IJ 

E8.3 

E8.3 

E8.  3 

E 13.  8 

E13.8 

E13.8 

13 

12 

Card  No.  3  EN  IK 

E8.3  12 

Card  No.  4  GAMMA  IL 

E13.8  12 

Card  No.  5  ELK  IM 
E8.3  12 

Card  No.  6 

E8.  3 
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Explanation  of  Symbols 


L 

NP 

KL 

II 

QT 


number  of  harmonics  to  be  calculated 

interval  of  printing  (+-  1  print  integral  for  every  llf  used,  +  2  for  every 
other  ,  etc) 

-  1  — *  print  table  of  SiPSI,  +  1  — >  suppress  fhis  table 
number  of  cards  No.  2  to  follow 

determines  limits  of  integration  (QT  =  -  «  will  cause  integration  from 
z/rf  ytoyi"t-yf2irm-  w/2  and  from  i/r  +  y  +  2  ir  m  -  tr/2  to  2  tr  m 

*¥+  y)  T 


EMU 

H 

AMBDA 

l 

EM 

m 

DEPSI 

interval  size  to  be  used  for  ~j>' 

DELPHI 

interval  size  to  be  used  for  in 

DELTA 

interval  size  to  be  used  for  in 

K 

number  of  values  of  to  be  used 

IJ 

number  of  values  of  n  to  follow 

EN 

n 

IK 

number  of  values  of  y  to  follow 

GAMMA 

r 

IL 

number  of  values  of  i.  to  follow 

m  + 


t 


*  r 

2  n  m  + 


t 


ELK  i 

IM  number  of  values  of  x  to  follow 

X  x 
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Output  Formot 


MU  = 


LAMBDA  = 


M  = 


N  = 


GAMMA  = 


L  = 


(optional  table) 

PSI  S1PSI 


S2PSI 


S3PSI 


S4PSI 


I 


S1PSI 

A-ZERO  = 

A  -  k  =  B  -  k  = 


w 


I 

I 


52PSI 

A-Zero  = 
A  -  k  = 


V 


S3PSI 

A-ZERO  = 
A  -  k  = 

S4PSI 

A-ZERO  = 
A  -  k  = 

l 


l 

i 
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Sarrple  Dota  for  Part  I,  B 


Card  No.  I 

Card  No.  2 

Card  No.  3 

Card  No.  4 

Card  No.  5 

Card  No.  6 

Card  No.  7 

Card  No.  8 

Card  No.  9 


10H  -1  H  -  .  15707963  E  *  1 

+  .  100  E  +  0  +  .  500  E  -  1  +  .  300  E  ►  1  f  .  13089969E  +  0  +■  .  13089969  E 
♦  0  f  .43633233  E  -  1  f  49  +  1 

+  .  300  E  +  1  +  2 

+  .  0000  0000  E  +  0  +  1 

+.  00  E  +  1  +  1 

+.  800  E  +  0 

+  .  20943950  E  +  1  +  1 

+  .  100  E  +  1  +  1 

+  .  800  E  +  0 
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PROGRAM  II.  ORIGINAL  STRAIGHT  LINE  APPROXIMATION  TO  WAKE 


This  program  attempts  to  define  the  trailing  wake,  "w",  by  means  of  the  very  much 
simplified  expressions  given  below(see  Ref.  1): 


w  = 


r 

TTT 


2  (  y  ~  7  )  cos  5 _ 

Z2  My  -  7  )2  cm2( 


except  for  ^  =  0  (1) 


^]/i  2  -  d2  sin  V'  (la) 


where  y  =  d  cosj 

♦  d=  M  [  5  +  f  -  <j>  ]  (lb) 

Z  =*15+/'-  ^  1  -aj/  -y)  (lc) 


+d  sin  7  sin  (  -  j'  ) 


tan  (  (£>  -  j/'  -  §  )  =  ^ 


-  u  cos  ^ 


(Id) 

de) 


i  f  H  sin^ 

Values  for  were  found  by  combining  equations  (lb)  and  (Id)  to  give  on  expression 

of  the  form  f  (  ^  )  =  0.0.  An  initial  guess  was  made  for  <f>  and  that  guess  improved 
by  a  Newton -Raphson  iteration  where  the  improved^equalled  the  original^minus  f(jzf)/ 

[f\  ^6  )  represents  the  derivative  of  f(  )  with  respect  to  <f>  ]. 


When  5  =  0 


Then  w  = 


r 

TTT 


1  -  sin  S 

(7-y)  cos  § 

whe.  o  ^  =  tan  ^  (  -b — c°s. 


since 


JL  +  M  sin  ^ 

<f>  =  f  .  d=  Z=0  y  =  Z 


This  program  also  performs  harmonic  analysis  of  the  function  as  summed  overthe  blades. 


The  function  4ttRw/  p  and  its  harmonic  coefficients  are  also  punched  out  in  a  format 
suitable  for  input  to  Program  III. 
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This  program  was  written  in  FORTRAN  II.  It  requires  approx imately  D5  minute 
on  an  IBM  7094  computer  to  handle  a  case  of  one^,  one  ^  ,  one  p,  and  yL  and 
4  C  at  intervals  of  15°  in  .  The  program  was  developed  at  the  Massachusetts 
Institute  of  Technology  Computation  Center  in  Cambridge,  Massachusetts. 

This  program  requires  two  subroutines: 

OUD1  -  to  compute  d  and^-  "Y* 

HANEW  -  to  perform  harmonic  analysis 

NOTE:  INPUT  AND  OUTPUT  are  taken  care  of  by 


WRITE  OUTPUT  TAPE  2 
READ  INPUT  TAPE  4  -- 


Therefore,  it  may  be  necessary  to  change  an  IOU  table  to  correspond. 

A  list  of  the  logical  to  physical  tape  correspondences  may  be  found  under  Program 
IA. 
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EXPLANATION  OF  SYMBOLS 


r 

R 

d 


1 

P 

X 

l 


a 

o 


s 

£ 


</> 


-  circulation 
=  rotor  radius 

=  horizontal  distance  travelled  by  rotor  hub 
=  vertical  distance  travelled  by  rotor  hub 
=  rotor  span  parameter  (equivalent  to  x  in  program  I) 

=  advance  ratio 
=  inflow  normal  to  rotor  disc 
=  rotor  azirrtith 
=  rotor  span  parameter 
=  coning  angl  e 
=  blade  spacing 

=  angle  between  vortex  line  and  a  perpendicular  to  the  blade 
=  wake  azirruth 
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INPUT  FORMAT 


Card  No.  1  II 
12 


Card  No.  2  EMU 
E8.3 


Card  No.  3  EL 
E8.3 


IK  NZ  NH  K 

12  12  13  13 


AMBDA 

E8.3 


ETA 

E8.3 


DEPSI  AO 

E 13.  8  E8.3 


IJ 

12 
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EXPLANATION  OF  SYMBOLS 


II  a  counter  to  be  set  equal  to  the  number  of  times  Card  No.  1  is  to  be  repeated 

IK  a  counte-  f  be  set  equal  to  the  number  of  times  Card  No.  2  is  to  be  repeated 

NZ  the  number  of  blades  to  be  considered 

NH  the  number  of  harmonics  to  be  considered  when  performing  harmonic  analysis 

K  the  nurrbei  of  values  of  ^at  which  the  function  is  to  be  computed,  must  be  odd 

DEPSl  the  interval  size  to  be  used  for  ^ 

AO  ~  a 

o 

EMU  =  M 

AMBDA  =  ^ 

IJ  =  counter  to  be  set  equal  to  the  number  of  times  Card  No.  3  is  to  be  repeated 

EL  =  / 

ETA  =  ^ 
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OUTPUT  FORMAT 


—  BLADES 


MU  = 


LAMBDA  = 


L  = 


ETA  = 


DELTA  PSI 


PSI  TOTAL  WAKE 

4  4 

HARMONIC  ANALYSIS 
A  -  zero  = 

A  -  i  =  B  -  i  = 


CONTRIBUTIONS  FROM  INDIVIDUAL  BLADES 

4 


Here  MU  LAMBDA  =  ?l  ,  L  =  J  ,  ETA  =  7  ,  PSI  =  ^  DELTA  PSI  =  interval 

;d  in  y,  . 


size  usee 


CONTRIBUTIONS  FROM  INDIVIDUAL  BLADES  -  These  blades  will  be  W 


4  TT  R 


A 


calculated  at  specific  blade  spacings  or  values 
of  .5  in  order  of  increasing 


TOTAL  WAKE  -  This  column  will  be  the  sum  of  W 


4  it  R 


7T 


over  all  blades  at  each 


P 


A  -  zero  ~  constant  term  resulting  from  harmo  c  analysis. 

A  -  i  coefficients  of  the  i^  cosine  te  '  suiting  from  harmonic  analysis. 


B  -  i 


coefficients  of  the  i^  sine  terms  resulting  from  harmonic  analysis. 


There  may  also  be  an  error  message. 

NO  CONVERGENCE  AT  ZETA  =  .  This  message  will  be  printed  out  if  there  is 

no  convergence  after  repeating  the  iteration  ten  times. 

In  this  case,  the  component  of  the  trailing  wake  under  consideration  is  set  to  zero. 
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PROGRAM  III  -  ROTOR  LOAD  ANALYSIS 


1.  Description  of  Procedure 

The  nondimensionalized  trailing  wake,  V  |  (  y  is  defined  over  a  range  of 

0  — ►  2tt  in  the  azimuth  angle  .  Values  of  ^  generally  vary  by  25°  or  15°.  For 

any  given  configuration,  "m"  values  of  and  "m  +-  1  H  values  of  £  ,  where  /  and 

^  are  rotor  span  parameters,  are  considered^o  that  there  are  "m^  t-  m"  sets  of 

V,  (^).  The  members  of  these  sets  are  designated  as  A. ^  ^ )  where  "  j "  refers  to  a 

specific  value  of  V  and  k  refers  to  a  specific  value  of  t  .  Harmonic  analysis  with 

2 

respect  to  y*  is  performed  on  each  of  these  sets,  so  that  there  are  also  "m  +  m"  series  of 

"2n  +  1"  harmonic  coefficients  P  and  Q  such  that  A.^  (  (y  )  ^  Pq  f  £] 

1  jk  <•/ 

(  P.  cos  i  +■  Q  .  sin  i  V ).  In  the  following  write  ■  up,  i  will  always  indicate 
'jk  1  jk 

a  harmonic,  j  a  specific  value  of  ,  and  k  a  specific  value  of  t  . 


A.  To  find  the  load  coefficients,  a  system  of  "m"  equations  is  set  up: 


o. 


b 

w 


1 1 

2  KvV(p0 ,-po  >1  ■ 


tan 


t  =  1 


jt  +  1 


and  the  ^  are  determined.  These  may  then  be  compared  with  the  "  A  "  or 

°i 

"  A  '"s  used  in  calculating  the  P's  and  Q's.  If  desired  an  iterative  process  may  be 

set  up,  replacing  the  "  A"  with  the  A  ,  until  reasonable  agreement  is  obtained  between 

i 

A  and  the  range  of  A 


*  P  *3  *  2nM 

v  (^/)  =  f  1  /  U  1  d  cos^-  Tcos  if-P)  +  H  sm^)  -p7sinH<^  /? 

(  i  fd  '  +1  f  Z  -2  'n  Jl  cos  ( f-f ')  +  2  /  d  cos^  -2  V  d  cos  j 


24 


A  A  and  a  may  also  be  obtained. 

lCj  is. 


m 

£ 

n 

t=  1 

m 

b 

W 

G, 

(™-WV. 

A  -  b 

A  ic.-  ^T 
I 


A  is.  * 

I 

B.  The  time  history  of  downwash  A  •  (  ^0  is  calculated  in  the  following  manner: 

m 

*j<W-  A0t)(Ait^)-Aittl(^))]  -H-ani 

An  alternate  method  [  A  •  (j^)  -  AQ  f  S  (  Ajc  cos  i^+  A.^  sin  i  | 

1  °i  .  'ci  ,si 

i  =  I 

was  rejected,  because  it  did  not  appear  to  give  a  very  accurate  reconstruction,  unless 
some  rather  high  harmonics  were  included. 

C.  The  distributed  airloads  Ln  s)  for  n  >  2  are  calculated  for  each^.  as  follows: 


Let  Y.  =  F.  A-  -G  .  A  • 

IC.  I  IC.  I  IS. 

I  1  I  I 


and  y.  =  F.  A  •  +  G.  A  • 
is.  I  is.  |  ic. 

I  I  I 


where  i  =  2,3,4,5,6,7and  Fand  G  are 
general ly,  but  not  necessarily,  constant 
over  i  and  j. 


Then  L. 

ic. 

I 


f  lYiti,., 


L. 


(c,s)j  for 


i  =  3, 4, 5,6 


“5*  1  ri  - 1,  c. 

•  i 


c. 


C  is  a  conversion  factor,  depending  upon  the  dimensions  of  the  test 
model  under  consideration. 
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D  The  I  cm  T  function  L.  (^)  is  computed  in  two  ports: 

(1)  L  L  (VO  ’  -c  (7  j  •  H  sin  If')  A|  (h 

(2)  Lq  (p)  =  '  c  [e.(^  .2  .  v7’7)  -7.  ho, 

i  *  2 

*  1 2  H0j7j  '  a]  (7j  “P  4)m2)I  sin^  +  [b  |  (7j2^p2  4)-p  ajl  i)  cos 
Mp  b]7  |  -  aQ)  sin  2  ^  *  (p  Oj  y.  -  -j-pV  )  COS  2/^ 

4  -4-  aj  sin  3  ^  -  -±j-  b,  cos  3  ^  } 

L:(^)=  L,  (yo  *  Lg  (^) 

i  i 

C  is  again  a  conversion  factor,  and  aQ,  Oj,  bj  are  geometric  parameters 
(blade  coning  angle  and  first  harmonic  flapping). 


L.  (y)  is  also  computed  with  the  zero~th,  first,  and  second  harmonics  eliminated. 

(The  third  harmonics  in  Lq  (^)  are  considered  negligible.) 

i 

L-  (V')  =  C^7-  •  H  sln^)  ^  .  (j^)  -7/1  Q  -  -V  Rl,. 

M7  ;  A  ]  f  7^  A  2s  '  cos^  "  (7j  A  lsi  *■  M  A0  -  A  2c  )  S'n  t* 

1  i  1  i  i 

-<7ihc.  '  '•  'V  ■‘j  *  3s.  >  cos  2  ^  "  (7j  ^2s.  ‘  *7  M.  '“^Xj  iin  2Zl 

I  1  I  II  I 

NOTE:  If  it  is  desired  to  do  a  case  wheie  the  number  of  I's  is  not  one  greater  than  the 
number  of  7  /  fben  a  different  procedure  must  be  used.  Suppose  Program  II  has  beenrun 
for  f  ive  values  of  7  and  hvo  values  of/. 

=  /M  P  “P  I  -ptani 

o.  /  O.i 

I  I  1 


j2 


X  •  =/>l  p  -  p. 

A  ic.  /  1..  i.0 

I  1 1  |2 


A  is.  Z5  1 . .  1 . 


Q.  -Q.  j 
‘i2 


and 


F  l-*-5 


A.  (j^)  lAji  0^)  ■  A  j2  ()^)  i  -  M  fan  i 
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Where  p  is  an  appropriate  conversion  factor.  A  rough  approximation  would  be 
b  times  the 

TIT 

These  "  ^  "  must  then  be  put  into  lie  appropriate  card  format  for  input  data 
to  Program  III  C.  If  a  small  computer  is  available  it  is  easy  to  write  a  simple  program 
to  do  this. 


average  .Q.~  ^  ) 
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2.  Description  of  Variables 


V|  (^)  =  non~dimensionalized  trailing  wake  for 
y*  =  azimuth  angle  for  rotor 
3  =  determined  by  blade  spacing 
M  =  number  of  wake  spirals 

l  l^l  ~  rotor  span  parameters  is  the  same  asiX  in  progran  I) 

d,  'Ll  -  horizontal  and  vertical  distances  travelled  by  rotor  hub 
(j  =  advance  ratio 
^  =  azimuth  angle  for  wake 

m  =  number  of  values  of  y  under  consideration 
n  =  number  of  harmonics  taken  in  harmonic  analysis  of  V|  0*) 
ic  or  is  =  indicates  an  ith  cosine  or  sine  harmonic  resulting  from  harmonic  analysis 


=  indicates  a  specific  value  of  y. 


F 

G 


=  indicates  a  specific  value  of  i-  J ^ 

=  lift  deficiency  function  -  from  data  obtained  from  shed  and  harmonic  trailing  wakes. 
=  lift  phasing  function  -  from  data  obtained  from  shed  and  harmonic  trailing  wakes. 


3.  Program  Mechanics 


This  program  is  written  in  FORTRAN  II.  It  requires  approximately  .  3  minutes 
to  handle  a  case  involving  five  values  of  y  on  the  IBM  7094.  One  subroutine  "CROUT" 
is  used  to  solve  the  system  of  simultaneous  equations.  A  maximum  of  7y,  0,  F  or  G,  a 
maximum  of  50  values  ^  ,  and  a  maximum  of  20  harmonics  may  be  used.  This  program 
was  developed  with  the  aid  of  the  facilities  at  the  Massachusetts  Institute  of  Technology 
Computation  Center. 
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INPUT  FORMAT 


Card  No.  1  DEPSI  BOER  TANMU 
E13.8  E8.3  E8.3 

EMU  CF  II 

E  8.3  E8.3  13 

M 

13 

KP 

13 

NH 

13 

NGP 

13 

I2P 

13 

Card  No.2  ETA  (1)  THETA  (1) 

E8.3  E8.3 

F(l)  G(l) 

E8.3  E8.3 

This  card  format  must  be  repeated  "m 

Card  No.  3 

"  times  in  succession 

A.  If  I2P  =  1 

1.  A(l,  J,  K)  5E14.7  1 

2.  BO,  J,  K)  6E12.  5  J 

group  ^ 

These  cards  are  to  be  grouped  by  ^  then  by  J  i.  e. 
en*er  i.  and  2.  for  followed  by  1  and  2  for  J  ^  etc-  ,  up 

,  hold  **1  constant; 

to  l.and  2.for  f  , 

m  +  1 

then  consider  ^  ^  and  repeat  the  process. 

B.  If  I2P  =  2 

In  this  case  the  ordering  will  be  the  same,  except  that  each  subgroup  of 
1  and  2  will  be  augmented  by  3  and  4 


1. 

A(l,  J,  K) 

5EI4.7 

2. 

B<l  ,  J,  K) 

6E12.5 

3. 

C(K) 

5E14.7 

4. 

D(K) 

6E12.5 

in  this  case 

1  and  3  are 

added  together  as  are  2  and  4 

Card  No.  4 

AZERO 

AONE  BONE 

3E8.3 
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Exploration  of  Voriobles 


DEPSI  =  interval  size  used  for  ^ 

B°ER=  £ 

TANMU  =  u  tan  i 
EMU  =  M 

CF  =  a  conversion  factor  depending  on  the  geometry  of  the  model  being  considered 
II  =  number  of  times  Card  No.  1  is  to  be  repeated 
M  =  number  of  values  of  ^ 

KP  =  number  of  points  of  ^ 

th 

NH  =  number  of  harmonics  (  o  term  +  all  cos  coefficients  +  all  sin  coefficients) 

NGP  =  number  of  card  No.  4  to  be  input 

I2P  =  +2  indicates  C  and  D  arrays  are  to  be  input,  +  1  indicates  only  A  and  B  are 

to  be  input 

ETA(I)  =  7. 

THETA  (I)  =  0. 

i 

F(D=  F. 

G(I)=G; 

A(l,  J,  K)  =  harmonic  coefficients  resulting  from  the  harmonic  analysis  of  B(l,  J,  K) 

B(l,  J,  K)  "KP"  values  of  the  trailing  wake  functiorv(lf  the  function  has  been  separated 
into  near  and  far  wakes,  this  will  be  the  far  wake  portion) 

*  C(K)  harmonic  coefficients  resulting  fromlhe  harmonic  analysis  of  D(K) 

D(K)  "KP"  values  of  the  near  wake  f.  nn  of  the  trailing  wake  to  be  added  on  to  B 

*  These  will  be  omitted  if  I2P  =  + 1 
AZERO  =  a 

o 

AONE  =  a  j  Blade  flapping  is  defined  as  aQ  -  a  |  cos  Ty  -  b  j  sin  ^  . 

BONE  =  bj 
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OUTPUT  FORMAT 


Page  1  LOAD  COEFFICIE NTS  FOR  ETA  =  ,  , 

HARMONIC  NO.  i  % )  (+\=1?2)  (  f  ( V  =  7W > 

V  i  +  1  il  i  v  Jz 

(These  will  follow  sequentially  so  that  the  first  sine  harmonic  will  be 
assigned  a  number  equal  to  that  of  the  last  cosine  harmonic  plus  one) 

Page  2  LAMBDA  (PSI)  FOR  ETA  = 
psi 

' l 


i 


*  .(P) 


Page  3  ETA  L3-C  L3-5  L4-C  L4-S  L5-C  L5-S  L6-C  LS-S 

J,  l  I  l  l  t  l  l  l 


Page  4  ETA  = 

PSI  L(l) 

l  I 


L(G)  L(l)  +  L(G)  L(G  +  |)  -  HARMONICS  EXTRACTED 

i  i  i 


Explanation  of  Symbols 
ETA  =  ^  PSI  = 

LOAD  COEFFICIENTS  ^  .  ,  A  is 

i  i  i 

LAMBDA  (PSI)  =  2  -(h 

Ln  -  C,  Ln~S  =  L.  ,  L. 

ic.  is. 


I 


UD  =  L,  (?) 


L(G)  =  Lg  (?) 


L(G  +  I)  =  HARMONICS  EXTRACTED  -  L.  (?) 
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APPENDIX  A 


Program  for  Vj 


Trailing  Wake 


CINE  N  S  ION  *  IPS  1 1  361  I  .  Y2PSI  1361  I  .V  JPSII  >611.  Y4PSI  <361  I  .Y8PSU  Hit  , 

1  TUP  1(3611.  »7PSI<3tll,»8PSil  36  1 1 , 2  APS  I  I  361 1 . 2BPSI  I  361  I ,PS I  (361  1  . 

2  CHl(36ll,SIGNA(36l),  AouAVI  201  ,  BBUA  Y  (  2C  I 

CONNCN  PS  I  .  T  IPS  I  ,  I . L ,x ,0E  PSI , AOUAV , BGUAY ,A2E»0, ▼ 2 P S I  . Y3PS I . 

1  r*PS  l  . Y5PSI , Y6PS I  . Y TPS I ,Y8P$  I 
1  REAO  INPUT  (APE  *  .  2  ,  l  ,AL .  NP,  I  1  ,NN  ,NN  ,GT 

2  FORMAT ( 612. E 1 3. 8 ) 

READ  1N»0T  TAPE  4.6.BC 

3  RfAC  INPUT  TAPE  *, 4, £nu, AN8CA,  Eh, CE P S I . C£ LPmI , CEL TA , A , I J 

a  roRPAt  1 3 e e . 3, 3E  i  3.e,  n,  12 i 

5  READ  INPUT  TAPE  4,6,A,I« 

6  f  CRhaT  I  f  e  .  3,  12  I 

T  REAL  Input  tape  a, 6. Ela.il 

8  KEAC  Input  tape  4.6.EN.IN 

9  K  f  A  C  INPUT  TAPE  *.6,*,  IN 

N  •  5 

10  -RITE  OUTPUT  TAPE  2 ,  1  1  .  E*U , AN80 A  ,  E N . E N , E L * , A , « , OE P S I , DE LPM I , PE L T A 

11  PCRPATIAH1PU  •  F  6 . 3  ,  1  3h  L  A  “  8  C  A  •  F 1,4, 9M  N  -  Fe.l, 

1  9H  N  ■  36.3,9m  L  •  F6.3.9M  A  •  F6.3.9H  X  •  f 6 .  3 / / 

2  I3H006LTA  PSI  •  F8.S.20H  CElTA  Pm  111  I  -  FB.S.20M  C  EL  T  A  Pm 

31(21  •  F  8 • 3 1 

CO  12  J  ■  1  .* 

CP  I  I J I  •  0.0 

5  I  UP  A  I J I  *0.0 

12  CONTINUE 

1)  REAO  INPUT  I  APE  A, 1*. GANNA 

14  FOftNAT  (£13.81 

N  A  N  -  1 

PSII1I  •  0.0 
CO  14  j  •  1  .  A 
ALPHA  •  0.0 
0"EGA  •  0.0 
SETA  -  0.0 
2  E  T  A  •  0.0 

PHI  •  PSIIJI  ♦  CANN* 

SPMl  •  P m |  .  6.2831852*14 

19  ASSIGN  2 3  TC  JA 

2C  0  •  EPU«(SPh|  -  PHI) 

l  ■  AN8CA«( SPhi -phi  I  •  I  1 .0*A*CCSF  I  PHI  I  1-0041 ELA-X  I 

RmO  •  ELA..24l»42*C«»24/»»2-2.0«>»£LA»CCSE  (  PSI  I  Jl -PHI  l*2.0*C« 

1  *  Ela»CCSF  IPh|  I  -  x  •  C  C  S  P  (PSII  J1  I  I 
IF  (RHCI  99,99,23 

23  TAU-IELA*(ELA«C«COSF(Ph|I-x»CCSEIPSIIJI-Ph|II«ENu"IEL«*SINF(PhII 
1  -X*S I NF ( PSI (J  III  I /(RhC** 1. S) 

uO  TO  JA, ( 25.  3  3.40,81 ,58,89  I 
25  ALPHA  •  ALPHA  ♦  TAU*C0SF (ENAPHl  I 

CNfOA  •  ONEGA  ♦  TAU«S  INF  I  EN«PH |  1 
PH |  •  nHl  A  CELPhI 

IFIGAnna  -.000011  29,29.43 

29  IF  ( Pm | - SPh | -t T  I  30,30,  50 

30  ASSIGN  33  TC  JA 
GC  TC  2C 

33  ALPHA  •  ALPHA  A  4.C  aTAUaCOSF  I  ENaPH  |  1 

CPC  G*  •  CNEG*  ♦  4.0aTAU»S INF  (  ENaPhI  I 
Phi  •  Ph I  a  CELPhI 

IF  I  PmI  a  0 . 5  aDE  L  Ph  1  - SPh I  I  35  ,  35  ,  1  9 
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35  IFIGAMMA-.OOOOII  37,37,  36 

36  ASSIGN  AO  TO  JK 
GO  TO  20 

37  1HPHI  40.5«OELPHI-SPH1-OT I  36,36,  19 

40  ALPHA  -  ALPHA  ♦  2 .0 • T AU*CO S F ( EN«PH I ) 

OMEGA  >  OMEGA  4  2.0«TAU» S INF ( EN»PHI ) 
PHI  •  PH!  ♦  DELPHI 
ASSIGN  33  TO  JK 
GO  TO  20 

43  IF  (PHI  -  SPHII  30, 30,45 

45  GO  TO  (71,69,67,65  I  ,M 

50  PHI  -  PH!  -  0ELPH1 

51  BETA  •  BETA  ♦  TAU*COSF( EN»PHl I 

ZETA  •  ZETA  ♦  T AU» S INF ( EN«PH I) 

PHI  •  PHI  ♦  DELTA 
IF  (PHI-  SPHII  52,52,63 
52  ASSIGN  55  TO  JK 

GO  TO  20 

55  BETA  «  BETA  ♦  4. 0« TAU» COSF ( EN«PHI ) 
ZETA  *  ZFTA  ♦  4.0»T AU»S l NF ( EN«PH l ) 

PHI  ■  PH!  ♦  OELTA 

IF  (PHI  ♦  0 • 5*DEL  T  A  -  SPHII  57,57,58 
57  ASSIGN  59  TO  JK 

GO  TO  20 

56  ASSIGN  51  TO  JK 

GO  TO  20 

59  BETA  -  BETA  *  2.0*T AU*COSF ( EN«PHI I 

ZETA  •  ZETA  ♦  2.Q4TAU4SINF(EN»PHII 
PHI  ■  Phi  ♦  DELTA 
ASSIGN  55  TO  JK 
GO  TO  20 

63  ZAPS  I ( J  I  ■  ( CELTA»BETAI /3.C 

ZRPSKJl  ■  (OELTA*ZETA)/3.0 
GO  TO  45 

65  Y1PSKJI  •  (DELPHI  .ALPHA)  '3.0 

Y2P9KJ)  •  (DELPHl«OMEGA  1/3.0 
CHKJI  ■  CHI  (  J)  ♦  VIPS!  I  J) 

SIGMA(J)  •  SIGMA(J)  4  Y2PSIIJ) 

GO  TO  73 

67  Y3PSKJI  •  (CELPHl»ALPHA 1/3.0 

Y4PSKJ)  •  (CELPHI»0MEGA  1/3.0 
CMIUI  ■  CHKJI  4  Y3PSI  ( J) 

SIGMA(J)  «  SIGMA(J)  4  Y4PSKJ) 

GO  TO  73 

69  Y5PSKJ)  -  (CELPHI«ALPHA  1/3.0 

Y6PSKJ)  ■  (OF.LPHI  *0MEGA)/3.0 
CHKJI  «  CHI(J)  4Y5PSI  (  J  I 
SIGMA(J)  •  SIGMA(J)  4  Y6PSIIJ) 

GO  TO  73 

71  Y7PSMJI  -  (OELPHI  4ALPHAI/3.0 

YPPSHJI  -  (OELPHI#OMEGA  I/3.C 
CHKJI  •  CHKJI  ♦  Y7PSI  ( J> 

S I GM A ( J  |  -  SIGMA(J)  4  Y8PSIIJI 

73  PS1 ( J4  1 )  -  PS  I ( J I  ♦  DEPSI 

74  CONTINUE 

IN  -  IN  -  l 
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!F(  IN)  77,77,  13 
77  1-3 

79  CALL  HARNAL  (PSl.CHlI 

WRITE  OUTPUT  TAPE  2,00 

80  FORMAT  (1H0////4H  CH  I /// / > 

CALL  PRINCE(L,AQUAY,BOUAY,AZERO) 

CALL  HARNAL  <  PS  I , SIGMA) 

WRITE  OUTPUT  TAPE  2,82 
82  FORMAT  (1H0////6H  SIGMA////) 

CALL  PRINC£(L,AQUAY,BQUAY,AZERO) 

CALL  HARNAL  (PSI.ZAPSU 
WRI TE  OUTPUT  TAPF  2, 04 
84  FORMAT  (IH0////6H  ZAPS!////) 

CALI.  PRINCE  (L»AQUAY,0QUAV,AZERO) 

CALL  HARNAL  ( P  S I , Z  BP  S I ) 

WRITE  OUTPUT  TAPE  2,85 
85  FORMAT  (1H0////6H  Zfa P S I / // / ) 

CALL  PRINCE!  L , AQUA Y , BQUAY , A ZE RC ) 

IF  (MM)  88,88,87 

07  CALL  TRIP  (CHI, SIGMA, PS1,ZAPSI,ZBPSI  ,NN,K) 

88  IF  ( KL 1  90,90,89 

09  CALL  PRICK  (NP) 

90  IM  «  IM  -  1 

IF  (IM)  92,92,9 
9?  IL  «  IL  -  1 

IF  (IL)  94,99,8 
94  IK  •  IK  -  I 

IF  (IK)  96,96,7 

96  IJ  »  IJ  -  1 

IT  (IJ)  98,98,5 
98  I  I  -  I  I  -  1 

IF  (ID  101,101,3 

99  WRITE  OUTPUT  TAPF  2 , 1 00 , PS  I ( ~ ) , PH  I 

IOC  FORMAT  (30HISINGULARITY  OCCURS  AT  PSI  OF  E16.0,8H  ,PHl  •  E16.8) 
GO  TO  90 

101  CALL  EXIT 

LND I  1 ,0,0,0, 0,0 ,0,0, 0,0, 0,0, 0,0,0 ) 
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bUPRCUllNE  HARNAL  C  PS  I  *  FPSI  I 

01  MENS  ION  PS  I ( 361) , FPSI (361 ) , AQUA Y ( 20 ) ,BOUAV( 20) 

COMMON  PS  I »FPS  i  • I • L »K*  OEPSI , AQUAY , BGUA Y , AZ ERO 

KN*K-3 

R-1.0 

GO  TO  (5,35,5)  ,1 

5  GRAL  «  FPSI ( 1 )  ♦  FPSI(K)  ♦ A . 0 -FP S I ( K- 1 ) 

DO  12  M»2  »  KN , 2 

12  GRAL  «  GRAL  ♦  A.C-FPSKM)  ♦ 2 . G -F PS 1( 1 ) 

AZERO  *  (OEPSI -GRAL )/9. 4247770 
(J«0. 

DO  25  J»1,L 
Q-G  *1.0 

GRAND  *  FPSI  (  1  )»COSF(Q»PSI  (  1  )  )  ♦FPSKK  )*COSF  ( Q-PS 1  (K)  ) 

1  ^.O-FPSI  (K-l  )*COSF(Q*P$I(K-l  )  ) 

00  2C  M«2,KN,2 

20  GRANO-GRANC  ♦ A  ,0-COSF ( Q »P S I ( M  )  )  «FPS I ( M ) 
l  ♦2.0*FPSI (M4l)»C0SF(G«PSI (M+l )  ) 

AGUAY(J)  -  ( ( 2.G-DEPSI )/9.4247778)-GRAND 
IF  (1-3)  25,22,22 
22  A GUAY (J) ■AQUAY(J)/2.0 
25  CONTINUE 

IF ( 1-3  )49, 34, 34 

34  AZERO  -AZER0/2.0 
R  «0.5 

35  P  -0.0 

no  43  J  - 1 ,  L 
P  «P  ♦l.C 

TGRAL  *  FPSI ( l  )-SINF(P*PSI(  111  ♦  F PSI ( K I • S INF ( P*PS  I  ( K )) 
1  ♦4.0-FPSI (K-l)«SINF(P»PSI(K-l)  ) 

DO  42  h  -2 , KN , 2 

42  TGRAL  -TGRAL  *4.0*FPS1(M)»SINF(P-PSI(M) ) 

1  *2.0 -FPSI (M41)«SINF(  P-PSI (*♦! ) ) 

OOUA Y ( J )  -(  (2.0 -0E PS  I )/9. 42 47778) -TGRAL «R 

43  CONTINUE 
49  RETURN 

END ( 1,0, C ,C, 0,0, 0,0, 0,0, 0,0, 0,0,0) 
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SUBROUTINE  PRICK(NP| 

OIMLNSION  YIP  SI ( 361) .Y2PSII361 ),Y3PSI (361),Y4PSI  I 36 1)  ,  Y5PS I ( 36 1 >  » 
1  Y6PS I ( 361 ) ,Y7PSI(361),  Y0PS I ( 36 1) , PS  I  ( 36 1 1 1 AQUAV ( 20  )  tBGUAY (20) 

COMMON  PS  1, VIPS  I,  I,l,K,OEPSI ,  AQUA Y , 8CUAY , A Z ERO , Y2P S I  ,Y3PSl, 

1  Y4PSI.Y5PSI , Y6PSI tY7PSl ,Y8PSI 
WRITE  OUTPUT  TAPE  2.3 

3  FORMAT ( 102H  IPS  I  Y1PSI  Y2PSI  Y3P 

IS1  Y4PSI  Y5PSI  Y6PS I /// ) 

Of-  9  J  »  l.K.NP 

WRITE  OUTPUT  TAPE  2 , 6  ,  PS  I  (  J  )  ,  Y  l  P  S I  (  J  »  ,  Y2P  S I  (  J  )  ,  Y  3P  S  I  (  J  )  , 

1  Y4PSI ( J), Y5PSI ( J) ,Y6PSI  (  J) 

FORMAT  (E13.5.6E16.5) 

CONTINUE 
RE  TURN 

ENC ( 1 »  Of  C  fOfOtCfOfOfOfOfOfOfOtOtO) 
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SUBROUTINE  T R IP ( CH I , S IGMA , PS  I , Z AP S I ,  BPSI,NN,K) 

DIMENSION  CHI (361 ) , SIGMA! 361 ) ,PSl (361) ,2  APS  I (3611 ,ZBPSI ( 361 ) 
WRITE:  OUTPUT  TAPE  2,1 

1  FOR  MAT ( 78H iCH I  SIGMA  PS! 

1  ZAPS  l  ZBPSI///I 

DO  9  J  ■  1 ,  K , NN 

WRITE  OUTPUT  TAPE  2 , 5 ,C H I ( J ) , S 1 GMA ( J ) , PS  I ( J ) , l APS  I ( J ) , Z BPS  I ( J  I 
5  FORMAT  (E13.StAE18.5l 

9  CONTINUE 

RETURN 

ENL( 1,0,0,0,0,0,0,0,0,0,0,0,0,0,01 
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9U8R0UT INE  PRINCF (L , AQUAY , 3QU AY  , AZERC J 
DIMENSION  AOU A Y ( 20  )  , BCUAY (20) 

WRITE  OUTPUT  TAPE  2, 5, AZERC 
5  FORMAT { IHO///IOHCA-ZERO  ■  E  1 0 . 8  > 

00  9  J  *  111 

WRITE  OUTPUT  TAPE  2 , 7 , J , AGUAY ( J ) , J , 8QU AY ( J  ) 

7  FORMAT  ( JHGA- 1 2  *  3H  ■  618.8, IOH  B-I?,TH  *  E18.8) 

*  CONTINUE 

RE  TURN 

E N C ( 1 ,0,C,0,0,0,0,0,0,0,C ,0  , C , 0 , 0 ) 
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APPENDIX  B 


Typical  Results  from  Program  of  Appendix  A 


*U  •  0.100  1**60*  •  0.0500  *  •  3.000  N  •  3.000  l  •  l.OCO  »  •  0.500  «  •  C.550 

CELT*  PS|  .  0.13090  cm*  Phllll  •  C.13C90  CELT*  PR  1  ( 2  I  •  0.0*363 


CM 


•  o.m«*MU  oi 


' 

• 

C . 23*3571 7( -uO 

6-  1  ■ 

O.I89?97*9E  01 

-c. 6*777776!  ci 

0-  7  • 

-C. 309*  75**E -00 

1  .  mo  70551  01 

ft-  1  • 

-0 . 1 8  8  7  9 3* IE -0  ? 

•» 

m 

-r.  37869507C  -00 

0-  *  • 

C  .  1  7 8  1  669  1 E -00 

-  •> 

■ 

C  .  7*906539f -0  1 

ft-  5  • 

-0.7*0  3 1 100E -00 

-  6 

• 

C. 773317501-01 

H  -  6  • 

0. 301 9979*E  -00 

-  7 

m 

-'  .  7  3  72  76061  -00 

6-  7  • 

-0 . 350C  5 79  IE  -00 

-  P 

• 

G.  1  1  *  1  56  3  1  1-00 

B-  8  • 

0. 71*3310 3C  -00 

-  9 

• 

-C. 337537761 -00 

B-  9  • 

•'.  770C  3673E-01 

-1C 

■ 

C.  3g  7*  *6*  31  -00 

P  -  1  C  • 

-C. 777*73?  5F-01 

SI'J** 


-7E«L  ■ 

C .*301 *069E -CO 

-  1  - 

-L'.**8t  1 1  3*f  -00 

ft-  1  • 

0.237*808*E 

01 

-  7  • 

0.6173  39081  00 

6-  2  • 

-C.26375015E 

01 

-  3  • 

-C.  1«90?7*H-O0 

6  -  3  • 

0.3*310731E 

01 

-  *  - 

0. 339566181-00 

6-  *  • 

-0. 166*978SE 

-00 
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A-  5  • 

-C. 206589626-00 

fi- 

5 

• 

0 .408243656-00 

A-  6  ■ 

-0.230542946-01 

B- 

6 

m 

-0.25071085E-00 

A-  7  a 

0.  385338786-00 

B- 

7 

m 

0. 122725826-00 

A-  8  • 

-0.416841666-00 

8- 

8 

m 

-0. l*012528E-0l 

A-  9  - 

0. 267635886-00 

e- 

9 

9 

-0.740385456 -01 

A-  10  • 

-0.213721046-01 

8- 

10 

m 

0. 16637603E-00 

ZAPSI 

A-ZERC  • 

-0.939 50 700E-03 

A-  1  - 

0. 2651 3033E-0  1 

8-  1  ■ 

O.59074956E-02 

A-  2  ■ 

-0.69089228E-01 

8-  2  ■ 

-0.24875500E-00 

A-  3  • 

O.61062555E  00 

B-  3  - 

0.154 3 1789E  01 

A-  4  a 

-0.21788670E-01 

0-  4  - 

0.572520366-02 

A-  5  • 

0.  1702  5375E-02 

8-  5  - 

0 . 278  75036E-02 

A-  6  ■ 

0. 42724407E-03 

B-  6  • 

-0.79660853E-04 

A-  7  a 

C. 169424646-04 

8-  7  a 

-O.53993321E-04 

A-  0  - 

-0.52219255E-05 

B-  8  • 

-0.53964967E-05 

A-  9  ■ 

-C .114500936-05 

B-  9  a 

0 . 262  1 1555E-06 

A-  1  C  - 

-0.220619966-06 

8-10  • 

0. 1785O340E-O6 

Z BPS  I 


A-ZERC  •  -O.2770V952E-O2 
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A-  1 

■ 

-C  .  576 1 290 4 E -02 

A-  2 

C 

G. 24882170E-00 

A-  3 

■ 

-G.15431759F  01 

A-  4 

■ 

-0. 57254104E-02 

A-  5 

a 

-0.27870707E-02 

A-  6 

a 

0*  80190378E-04 

A-  7 

a 

C. 54488785E-04 

A-  8 

a 

0 •  58944438  E-05 

A-  9 

a 

C . 2 146  1810E-Q6 

A-10 

a 

C.  3547  3040E-06 

B-  1 

S 

0.271 1821QE-01 

B-  2 

a 

-0.69089957E-01 

8-  3 

a 

0.61861937E  00 

B-  4 

■ 

-0.21786939t-01 

0-  5 

a 

0. 17027783E-0? 

e-  6 

a 

0.42748532E-03 

e-  7 

a 

0.1711481  IE “04 

B-  8 

a 

-0 • 5000  581  IE -0  5 

B-  9 

a 

-0.955 12301E-06 

B-  10 

m 

-0 . 1 30  3  85 1 5E “07 
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ChJ 

-0. 14443E 

01 

-0. 1 1941 E 

01 

-0. 76O10E 

00 

-0. 12395F 

-00 

0  •  7 1  14  1  E 

00 

0. 16927E 

01 

0 . 2  7  8  30  E 

01 

C.40208E 

01 

C.54198E 

01 

0.69236E 

01 

0. 04441E 

01 

0.97849E 

01 

0. 10868E 

02 

0. 1 1421E 

02 

C.  1  1685E 

02 

0. 10882C 

02 

O.90817E 

01 

0  85090E 

01 

0. 71* 90E 

01 

0. 52338E 

01 

0. 14907E 

01 

-0. 38389E 

01 

-0. 77273C 

01 

-D.53550E 

01 

-0. 55942E 

01 

-0.46113E 

01 

-G.64174E 

01 

-0. 35379E 

01 

-0. 14294E 

01 

C.48191E 

-00 

0.21510E 

01 

0.48529E 

01 

O.65OO0E 

01 

0.71 178E 

01 

0.75489E 

01 

0. 75093E 

01 

0. 71530E 

01 

0.631  136 

01 

C. 52503E 

01 

C.40534E 

01 

0. 28464E 

01 

C. 16618F 

01 

C.59144E 

00 

-0.26154E- 

-00 

-O.05951E 

00 

-0. 12390E 

01 

-0. 14528E 

01 

-0. 15240E 

01 

-C. 14443E 

01 

S  I  0*-i  A 


0.89760E  00 
0.18391E  01 
G .  250041  01 
C.27826E  01 
0.26691E  01 
C.22028E  01 
C.14463E  01 
0.48433E-00 
-0.53830E  00 
■0.1  340  6 £  01 
-C.17467E  01 
-0.17127E  01 
-0.96056E  00 
-0.68259E-G1 
0. 1631 7E  01 
0.34286E  01 
C.42460E  01 
0  4  8  8  8  9  E  01 
G.48419E  01 
0.6569 7E  01 
0.85724c  01 
0.10007E  02 
0.74.171E  01 
0.37151E  01 
C.  16  7  8  3  E  01 
-0 . 20  34  5E -00 
-0. 3229  IE  01 
-0  .‘•920  7E  01 
-0.6573CE  01 
-0.61685E  01 
-0.71092E  01 
-O.50772C  01 
-0.33759E  01 
-0.18400E  01 
-0.10496E  01 
-0.24684E-00 
0.44548E-00 
0.10630E  01 
0.11142E  01 
0.81138E  00 
0. 27679E-00 
-0 . 4  36 1 8E -00 
-0.11775E  01 
-0.17743E  01 
-0 . 20126E  01 
-0. 1 7726E  01 
-0.11042E  01 
-0. 15184E-00 
0.89760E  00 


P  S  I 


c. 

0. 13090E-00 
0.26 1B0E-C0 
0.3927CE-00 
0.52360L  00 
C.65450E  CC 
0.76540E  00 
0 . 9 1 6  30C  00 
0.1C472E  01 
C.11781E  01 
0.13090E  Cl 
0.14399E  01 
0.157C6E  01 
0.17C17E  Cl 
O.10326E  Cl 
C.19635F  01 
0.20944E  Cl 
0.22253E  01 
0.23562E  Cl 
C.24071E  01 
0.26180E  Cl 
C.274B9E  01 
C.28798E  Cl 
0.3C107E  01 
0.3M16E  01 
0.32725E  01 
0.34G34E  01 
0.35343F  Cl 
0.36652E  Cl 
0 . 3796  IE  Cl 
0.39270E  01 
0.40579C  01 
0.418B8E  01 
0.43197E  Cl 
0.44506E  Cl 
G.45815E  Cl 
0.47124C  01 
0.48433E  01 
C.49742E  01 
0.51051E  01 
0.52360E  01 
C.53669F  01 
0.54976L  Cl 
0.56287E  01 
0.57596E  01 
0.58905E  01 
0.6C214E  Cl 
0.61 52 3E  01 
0.62832E  01 


l  APSI 


0.55547E  OC 
0.  1C44  3 f  0 l 
0.  1  36  7  7E  01 
..U7G0E  01 
0.  1  3335E  01 
C.97399E  OJ 
C.44572E-00 
-0.1  7012E-0G 
-0.77732E  CC 
-0. 1 2795E  01 
-C.  15951E  01 
-0.  1  6  705F  01 
-0.14885E  01 
-0. 10714E  01 
-O.47042E-GC 
0. 20  ^  5  1  E -CC 
0.87248E  CC 
0. 14273E  01 
0.  1  70286  01 
0.  108  3  56  0  1 
0. 1 7115E  01 
0 . 1 290  1  E  01 
G.67972E  OC 
-0. 30368E-O 1 
-0.73B25E  OC 
-0.134C0E  01 
-0.  1  7500 E  01 
-0. 191 1 4 E  01 
-0.18045E  01 
-0.14492E  01 
-0.90139E  OC 
-C. 24385E-0C 
0.42639C-0C 
0. 10126E  01 
0.14323E  01 
0.16292E  01 
0.l5804t  01 
0.1 299  3E  01 
0.83325E  OC 
0. 25629E-0C 
-0. 341 84E-0C 
-0. B6A86E  OC 
-0.12438E  01 
-0.14093E  01 
-0.13406E  01 
-0.10496E  01 
-0.58404L  00 
-0. 18457E-01 
0.555461  OC 


ZHESI 


-0.131136  01 
-C.98P09E  00 
-0 .49453E-0C 
0.913976-0!. 
0.67M1E  OC 
0.116516  01 
0.14779E  01 
0.15609F  01 
0.13S57E  01 
0.10C32E  01 
C.43979E-C0 
-0.2 1CU6E-00 
-0.84722E  CO 
-0.13728E  01 
-0.17C36E  01 
-0.178516  01 
-0.16C02E  01 
-C.11721E  01 
-C.56127E  00 
0.14352E-00 
0.83875E  00 
0.14218E  01 
0.18C64E  0! 
0.19361E  01 
0.17921E  01 
0.13965E  01 
0 . 80624E  00 
0.11413E-00 
-0.58468E  00 
-0.11879E  01 
-C.16111E  01 
-0.17975E  01 
-0.17262E  01 
-0.14140E  01 
-0.912 54 E  00 
-0.29957E-00 
0.33241E -00 
0.88990E  00 
0.12919E  01 
0. 14813 E  01 
0.14336E  01 
0.11600F  01 
0.70621E  OC 
0. 14639E-00 
-C.43253E-00 
-0.93490E  00 
-C.120O0E  01 
-0.14138E  01 
-C.13113E  01 
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p  S I 


J. 

U.  I  >C»0£-PC 
C. 761801-00 
0.  192707-00 
0.52)607  00 
0.656507  00 
0.  7  8  5607  00 
0. <716107  00 
0.106777  01 
0.117817  01 
0.  1  >0907  01 
0.161997  01 
C. 15  7087  01 
0.17C  1  77  0  1 
0.181767  01 
0.196157  01 
0.709667  01 
C. 777517  01 
0.715677  01 
0.768717  01 
0.761801  01 
0.776897  01 
0.787987  01 
0.101077  01 
C.  116167  0  1 
C.  177757  01 
0.  160167  01 
0.151617  01 
0 . 166577  0  1 
0.179617  01 
0.  197707  01 
0.6C5797  01 
0.618887  01 
C. 611977  01 
0.665067  01 
0.658157  01 
C. 6  7  1767  Cl 
C. 686117  01 
0.697677  01 
( .510517  01 
0.571607  0  1 
0. 5 >6697  01 
C.  569787  01 
0.567877  01 
0.  5  75967  0  1 
0.589057  01 
C. 607167  01 
0.615717  01 
0.678  177  0  1 


V  IPS  I 


0.669797 -01 
0.187677-00 
C. 167697 -00 
0.516717  0C 
0.68  7677  0C 
0.815167  OC 
0.107107  01 
0. 1  11767  01 
0. 1  76  1  77  01 
0.767767  01 
0.17C717  01 
0.185097  01 
0.619867  01 
0.609077  01 
0. 136107  01 
0.186867  01 
0.959677  00 
0.119197  01 
0.171707  01 
0.501177  01 
0.55509L  01 
0.600707  01 
0.666597  01 
0.676887  01 
0.576387  01 
0.557807  01 
0.577567  01 
0.667907  01 
0.615757  01 
0.175777  01 
0.297777  01 
0.209767  01 
0.160857  01 
0.6  711  77  OC 
0. 696567-00 
0.785777  00 
0.178667  01 
C. 153757  01 
0.156667  01 
0.160787  01 
0.118167  01 
0.956)67  00 
0.76)927  00 
0.567677  00 
0. 167707-00 
0. 1 68897-00 
0.617127-01 
0.5 13017-02 
0.669797 -01 


T2PSI 


0.268657-00 
0.656077-00 
0.569  127  00 
0.650787  00 
0.751107  00 
0.885987  00 
0.960607  00 
0.671  167  00 
-C.  121977  -00 
-0.106627  01 
-0.129667  Cl 
-0.570277  00 
C. 852617  00 
0.260  117  01 
0.625607  01 
0.627267  01 
0.167027  01 
-0.182  607  01 
-0.691107  01 
-0.508517  01 
-0.196907  01 
-0.299667  01 
-0.116557  01 
-0.517177  00 
0.671627 -00 
0.1  10827  01 
0.226957  01 
0.101217  01 
0.157167  01 
0.610517  01 
0.686997  01 
0.666667  01 
0.621907  01 
0.298617  01 
0.850627  00 
-0.972117  00 
-0.192077  01 
-0.192627  01 
-0.151697  01 
-0.981187  00 
-0.566157  00 
0.311  867-00 
-0.266617-00 
-C. 107577-00 
-0.156587-00 
-0.119667-00 
-0.180057-00 
0.161167-01 
0.268857 -00 


7  3P51 


C.  52*107  -01 
-0.  nil  17-01 
-C. 160797-00 
-C. 2115)7-00 
-0. 160087  -00 
0.5  797  17-01 
0.652627-00 
0.986127  00 
C. 161957  01 
C. 22*907  01 
C. 281117  01 
C.  1085*7  01 
C.  101667  01 
0.291627  01 
C.16C117  0! 
0.617057  01 
0. 681007  01 
0.671917  01 
C. 602711  01 
0.287857  01 
0.11*667  01 
-0.271*87-00 
-0.1  7*167  01 
-0.293707  01 
-0.1**717  01 
-C.  18  3787  01 
-C.16C617  01 
-0.266267  01 
-C.  1  12867  01 
0.28**27-0*' 
C. 159897  01 
0.266577  01 
C. 1)5217  01 
0.158967  01 
C. 166927  01 
C. 296807  01 
C. 211697  01 
C.  160687  01 
0.117767  01 
0.7090  17  00 
0.66 3097  -  00 
0.1961  77-00 
0.160517  -01 
-0.5*6917-01 
-0. 501967-01 
0.  J61187-02 
C.  616167 -01 
0.8621  67-01 
0.52*117-01 


7  6  P  S  I 


0.6728*7-00 
0.871*77  00 
C.  122  167  0  1 
C.  167117  01 
0.1526*7  Cl 
C.  121567  01 
0.56  1077  OC 
-C. 100827-00 
-0.813207  00 
-0. *97167  00 
-0.121867-00 
0.619267-00 
0. 166117 -CO 
-0.110087  01 
-0.26**57  01 
-0.308617  01 
-0.215011  01 
-0.611157-00 
C.  120127  01 
0.269927  01 
C.  112217  0  1 
C.15C807  01 
0.10*957  01 
0.21*717  01 
C. 936617  00 
-0.6*9*71-00 
-C.  1*6977  0  1 
-0.289977  01 
-C.  165097  0  1 
-0.361777  01 
-0.2915*7  01 
-0.19*867  01 
-0.712917  00 
C. *01517  00 
C.  1  79977  Cl 
0.272607  01 
C.27697t  01 
C.  222997  0  1 
l.  1  18167  01 
0.1*6717-00 
-0.6*9197-00 
C.  809)51  OC 
-0.901797  00 
-C. *7)517  00 
-0.7*1217  00 
-0.56  )677  00 
-0.2*07)7  -OC 
0.71*227  -0  . 
C. 672*67  -I  0 


V  5PS  I 


-0.156197  Cl 
-0.1)65*7  01 
-0.9*1*87  CO 
-0.66*707-00 
0. 186067-00 
0.799*07  00 
0.1)0767  01 
0.172197  01 
0.20  1**7  01 
0.22)207  01 
0.262377  01 
0.2*6*57  01 
0.16)1)7  01 
0.661657  01 
0.692097  01 
0.685507  01 
0.609227  01 
0.2)7*07  01 
-0.900767-01 
-0.2*7587  01 
-0.562697  01 
-0.956927  01 
-0.126327  02 
-0.81**77  01 
-0.7*8787  0  1 
-0.6)0157  01 
-0.128*67  01 
-0.577657  01 
-0.6658)7  Cl 
-0.15*617  01 
-0.262017  01 
0.965827-01 
0.176027  01 
0.2855)7  01 
0.160)27  01 
0.177*17  01 
0.155)57  01 
0.108907  01 
0.2506)7  01 
0.185*67  01 
0.120117  01 
C. 511067  00 
-C.  1*85)7-00 
-0.750)27  00 
-0.115707  01 
-C. 1*1227  01 
-C.  155*07  01 
-0.1*0297  01 
-0.156197  01 


y**5  I 


0 . 1  55  897  - OC 
0.512*07  CO 
0.707*17  00 
0.650*07  CO 
C.19I  1*7-00 
C.*l  1167-01 
-0. 171907-01 
0.1 1 1*17-00 
C.  )9**7|  -OC 
0.602967  CO 
-0.12*6*7-00 
-0.15*167  01 
-C. 2  1  5  7  17  0  1 
-0.1559*7  Cl 
0.62  1  677  -0  1 
0.226207  01 
0.672*27  Cl 
0.712*17  01 
0.857167  01 
0.915557  01 
0.919917  01 
0.969677  01 
0.*19)lf  01 
C.2C652F  01 
0.270507-00 
-0.102207  01 
-0.1*2*97  01 
-0.506127  01 
-0.5*95*7  Cl 
-0.70)591  Cl 
-0.906)57  Cl 
-0. *57527  01 
-0. *90207  Cl 
-0.5629*7  01 
-0.1*9977  01 
-0.1*9*77  01 
-0. 18  1677-00 
0.7571*7  00 
0.166  757  0  1 
0.1629*7  01 
0.11)017  01 
0.7050)7  00 
-0.9271)7-02 
-0.59)267  00 
-0.89*777  00 
-0.8*96)1  00 

-0.66)607  00 
-0.2*1977  -00 
0. 155*97 -00 
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Hu  •  0.200  l»PBO*  •  C.0250  x  •  1.000  N  .  0.  I  •  C.'jOC  A  •  c.  *  •  0.850 

L'fL  T*  P8|  •  0.26180  Cllt*  PHII1I  •  0.26180  Cl  l  I A  Pp  I  I  2  I  •  0.0*16) 


O  I 


*-2fKC  •  0.2*896»50t  01 


I 

0.T6TJJM Jf 

01 

B-  1 

-0.592C90*« 

01 

2 

G. 959716921 

01 

8-  2 

-0 . *80  3*87  If 

01 

0.  92279*891 

01 

B-  1 

-0 . 2**8286  If 

01 

4 

G.  7001  8*071 

01 

8-  * 

-  0  •  6  3*096  l  2€ 

00 

} 

0.  180966681 

01 

B-  5 

0 .598*691*1 

00 

6 

i,.  58668*7*1 

OG 

6-  6 

0.1 16*J<»5  3E 

01 

1 

.  207511811 

01 

B-  7 

0. 1  1  170828' 

01 

fi 

-G. 17C507551 

01 

p  -  a 

0.6 1*6285  If 

00 

9 

-C. *22111621 

01 

8-  9 

-0. 176C2880E 

-00 

10 

-i  .  372632821 

01 

P  - 1 0 

-0. 109  1 1  5 1  *1 

01 

ill."* 


-2f80  • 

-  I  •  0. 

-  2  ■  -0. 

-  1  •  -0. 

-  *  •  -0. 


8-  1  •  -0. 

B-  2  •  -0- 

e-  i  •  -o. 

B  -  *  •  0  . 
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A-  5 


-0. 

i-  6  ■  0. 

i-  7  ■  0. 

-  0  ■  0. 

-  9  •  0. 

-10  •  -0. 


ZAP  5  I 


WERC  • 

-0. 1AA11066C 

A-  1 

■ 

-0. 50505327C  00 

A-  2 

■ 

-0.  1A080133E-00 

A-  1 

■ 

-0. 35666509E-01 

A-  A 

■ 

-0. 73993352E-02 

A-  5 

■ 

-0. 10136589E-02 

A-  6 

■ 

0. 55ABAA71E-0A 

A-  7 

a 

0. AA21 7125^-03 

A-  0 

m 

O.25120OA7E-O2 

A  J  9 

a 

0. 119O3009E-O1 

A-  10 

■ 

0. A69381C9E-01 

ZOPSl 


B- 

5 

a 

0 

B- 

6 

a 

0 

0- 

7 

a 

0 

B- 

8 

a 

-0 

B- 

9 

a 

-0 

B- 

10 

a 

-0 

B-  1  • 

-0.5651965AE  00 

B-  2  ■ 

-0.58095719E-01 

B-  3  • 

O.AOOC8912E-02 

B-  A  ■ 

O.5A272130E-O2 

B-  5  ■ 

0.230075B2E-02 

B-  6  ■ 

O.09713659E-O3 

B-  7  - 

0.902 151A5E -03 

a-  e  • 

0. 18273765E-02 

e-  9  ■ 

0.133 1A038E-02 

e-io  • 

-0.19368302E-01 

A-ZFRC  •  -0. 
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A-  l 


-0. 


8-  1  - 


0. 
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CHI  SIGMA 

0.35169E  02  0. 

0.27842E  02  C. 

-0.47245E  01  0. 

-0.21978E  02  0. 

-0.13122E  02  0. 

-0. 82261 E  01  0. 

-0.56000E  01  0. 

-0.39901E  01  C. 

-0.29033E  01  0. 

-0.21134E  01  0. 

-G.15026E  01  0. 

-0.1O037E  01  0. 

-0.87500E  00  0. 

-0.1 8872E-00  0. 

0.17  5 12E-00  0. 

0.53181E  00  0. 

0.89274E  00  0. 

0 . 1 26 29 E  01  0. 

0.16290E  01  0. 

0.19150E  01  0. 

0.18124E  01  0. 

0.95603E  OO  0. 

0.21277E  02  0. 

C.33543E  02  0. 

0.35169E  02  0. 


PS  1 

2  APS  I 

2  BPS  I 

0. 

-0.21508E 

01 

-0. 

O.2610OE- 

■00 

-0.22650E 

01 

-0. 

0.52360C 

00 

-0.22878E 

01 

-0. 

0.78540E 

00 

-0.22422E 

01 

-0. 

0.1C472C 

Cl 

-0.21506E 

01 

-0. 

0.13090E 

01 

-0.203C7E 

Cl 

-0. 

0. 15708E 

01 

-0.18949E 

O’. 

-0. 

0. 18326E 

01 

-0.17515E 

01 

-0. 

0.20944E 

01 

-0.16060E 

01 

-0. 

0.23562E 

01 

-0. 14624E 

01 

-0. 

0.26180E 

01 

-0. 13233E 

01 

-c. 

0.28798E 

01 

-0. 11910E 

01 

-0. 

0.31416E 

01 

-0. 10677E 

01 

-0. 

0.34034E 

01 

-0.956C2E 

00 

-0. 

0.366S2E 

01 

-0.85945E 

00 

-0. 

0.39270E 

01 

-0.7828 7 t 

00 

-0. 

0.M888E 

01 

-0.73343E 

00 

-0. 

0.44506E 

01 

-0.72136E 

00 

-0. 

0.47124E 

01 

-0.76070E 

00 

-0. 

0.49742E 

01 

-0.86B27E 

00 

-0. 

0.52360E 

01 

-0. 1058  IE 

01 

-0. 

0.54978E 

01 

-0.13276E 

01 

-0. 

0.57596E 

01 

-C.16415E 

01 

-0. 

0.60214E 

01 

-0. 19354E 

01 

-0. 

0.62832E 

01 

-0.21508E 

01 

-0. 
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>51 

T18SI 

7285  1 

V38SI 

Y4P5I 

Y  5P$  1 

Y6PSI 

0. 

0,578116  01 

0. 

0.237966 

02 

C. 

0.609176 

01 

0. 

0.261806 

-00 

0. 803  7  IE  01 

0. 

0.118236 

02 

c. 

C. 796206 

01 

0. 

0.523606 

00 

0 . 4 1 082E  01 

0. 

-0.189736 

02 

c. 

C.  101416 

02 

0. 

0.785406 

00 

-0.4  3  700E  01 

0. 

-C.  103396 

02 

0. 

-0.731936 

Cl 

0. 

0.104726 

01 

-0.22420E  €1 

0. 

-0.642866 

01 

G. 

-0.445146 

Cl 

0. 

0. 1 30906 

01 

-0. 1 2 7  1  4E  01 

0. 

-0.441776 

01 

c. 

-0.253706 

01 

0. 

0.  157086 

01 

-0.77806E  00 

0. 

-0.322796 

01 

c. 

-0. 15940E 

01 

0. 

0.183266 

01 

-0.485 J7E-C0 

0. 

-0.245456 

01 

0. 

-0. 10502E 

01 

0. 

0.209446 

01 

-0.291266-00 

c. 

-0.191666 

01 

0. 

-0.695466 

CO 

c. 

0.7 35626 

01 

-0. 1 5C21E-C0 

0. 

-0. 152196 

01 

c. 

-C. 44  1  346 

-00 

c. 

0.261806 

01 

-0. 38865E-01 

0. 

-0.121886 

01 

c. 

-0.244976 

-00 

0. 

0 . 2  8  7  98  6 

01 

0.56343E -01 

0. 

-0.976876 

00 

c. 

-0.8311 06 

-01 

0. 

0.  314166 

01 

0. 14438E-C0 

0. 

-0.7771 3E 

00 

0. 

0.577526 

-01 

0. 

0.  340346 

01 

0.23212E-00 

0. 

-0.60  75  36 

00 

c. 

C.  186696 

-oo 

0. 

3. 366526 

01 

0.325836-00 

0. 

-C  .460546 

-00 

c . 

C.  309836 

-00 

0. 

0. 392706 

01 

0. 432226-00 

0. 

-C. 337316- 

-00 

c. 

0.431906 

-"0 

0. 

0.4  1  8866 

01 

C.55941E  CO 

0. 

-0.223236 

-00 

0. 

0.556566 

00 

0. 

0.44  5061 

01 

0.71 8056  OC 

0. 

-0.141366- 

-00 

c. 

0.686226 

00 

0. 

0.471 246 

91 

0.922396  00 

0. 

-0.113566 

-00 

c . 

0.820176 

00 

0. 

0.497426 

01 

0.119086  01 

0. 

-0.224146- 

-00 

c. 

0.948396 

00 

0. 

0.523606 

01 

0.154486  01 

0. 

-C.  766006 

OC 

c. 

0. 103366 

01 

0. 

0.549786 

01 

0.2C1156  01 

0. 

-0.707706 

01 

c. 

0.102156 

01 

c. 

0.575966 

01 

0.265616  01 

0. 

0.166416 

02 

0. 

0.  197946 

01 

0. 

C. 602141 

01 

0.364976  01 

0. 

0.23979E 

02 

c. 

0.591446 

01 

0. 

0.628326 

01 

0.528136  01 

0. 

C.23796E 

02 

0. 

0.609176 

Cl 

0. 
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APPENDIX  C 


Program  for  ^ 2  ~  Wake 


01  MENS  ION  PSI ( 361  I, S1PSI ( 361 ) ,S2PSl ( 361 > ,S3PSI (361 )  ,S4PSI  (361  1  , 
1  AQUAY 1231  , BOU A Y ( 20  ) 

COMMON  PSl.SlPSI.l.L.K.DEPSI , AQUA V ,BQU AY , A? ER 0 
READ  INPUT  TAPE  4 , 3  ,  L,  NP , KL  ,  I  I  ,  QT 

3  rORMAT  (412, E13.fi) 

4  READ  INPUT  TAPE  4 , 5 , EMU , A MBOA, EM , DEP S I , OE IPH I , DEL T A , K , I J 

5  FORMAT  ( 3F8.3, 3E13.8, 13, 121 

6  RE  AO  INPUT  TAPE  4,7,EN,|K 

7  FORMAT  I E8. 3 » 121 

0  RE  AO  INPUT  TAPE  4, 9, GAMMA,  IL 

9  FORMAT  (E13. 8,121 

10  RE  AO  INPUT  TAPE  4,7, ELK, 1M 

11  READ  INPUT  TAPE  4, 12, X 

12  FORMAT  (E8.3I 

PSI (11  -0.0 

WRITE  OUTPUT  TAPE  2 , 15, EMu.AMBDA , EM, EN, GAMMA , ELK, X 
IS  FORMAT ( 6H 1  Mu  -  F6.3.12M  LAMBDA  ■  F6.3.7H  M  -  F6.3.7M  N  • 
l  F6.3.11M  GAMMA  »  F0.5.7H  L  -  F6.3,7H  X  ■  F6.3/////1 
00  43  J  *  1 , K 
ALPHA  *  O.j 
BETA  *  0.0 
ZETA  -  0.0 
ETA  «  0.0 

SPHI  «  6.2831852*EM*PSI  (  J I  ♦  GAMMA 
EPHI  ■  SPHI  ♦  QT 
PHI  -  PSKJ1  ♦  GAMMA 

17  ASSIGN  22  TO  JJ 

18  D  •  E  MU* ( SPHI  -  PHII 
Z  .  A MB DA • ( SPH I-PH I  1 

P!  *  Z**2MX*S  INF!  PSI  (  JI-PHP  *0*  S  I NF  1  PH  I  )  I  ••  2 
IF  (PI)  90,90,20 

20  RHO  -  x**2*Z**2*0**2-2.0*0*X*COSF(PSI( Jl) 

SIGMA  ■  RH0*ELK**2-2.0*X*ELK*C0SFIPSI ( JI-PHI ) *2.0*D*ELK*C0SF(PH| ) 
THETA  «  0»COSF ( PH  I  I  -X*COSF ( PS H J )-PH| ) 

TAU  ■ ( ( X  *S INF ( PSI IJI-PHI M0-SINFIPHI1 )/PI ) • I ( ELK* THE T A  1 /SQRTF 
1  (SIGMA)  -  THET A/SQRTF(RHO|  I 

GO  TO  JJ  ,  (22,25,30,35,40,45) 

22  ALPHA  *  ALPHA  ♦  T AU *C0 SF ( E N*PH I) 

BETA  -  BETA  ♦  T AU* S I NF ( E N*PH I  ) 

PHI  «  PHI  ♦  DELPHI 
ASSIGN  25  TO  JJ 
IF  (PH1-EPMI)  18,23,23 

23  PHI  -  PHI  -  DELPHI 

24  ASSIGN  35  TO  JJ 
GO  TO  18 

25  ALPHA  «  ALPHA  ♦  4.0*TAU*C0SF( EN*PH|) 

BETA  *  BETA  ♦  4.0*TAU*S INF( CN*PHI ) 

PHI  *  PHI  ♦  DELPHI 

I F ( E  PHI  -(PHI  ♦0.5*UELPHII)  1*, 17,27 
27  ASSIGN  30  TO  JJ 

GO  TO  18 

3u  ALPHA  ■  ALPHA  ♦  2 .0*TAU*COSF(EN*PHI  ) 

BETA  -  BETA  ♦  2.0*TAU*SINF(EN*PH| ) 

PHI  -  PHI  ♦  DELPHI 
ASSIGN  25  TO  JJ 
GO  TO  18 
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35  ZETA  -  ZETA*TAU«SINF(EN*(PSI(  Jl-PHI) )»COSF(EN*PSI I  J)  ) 

ETA  -  ETA  ♦  TAU«S IMF ( EM* ( PSI I J )-PHl ))*SINF(EN*PS((J)) 

PHI  ■  PHI  ♦  DELTA 

ASSIGN  AO  TO  JJ 

IF  (PHI  -SPHI)  18* SO*  50 

40  ZETA  -  ZETA  f 4.0*TAU»S INF ( EN* (PSI  (  J  t-PHI  )  )  *COSF  ( EN«PS1  ( J)  ) 

ETA  ■  ETA  ♦  4.0«TAU«SINHEN*(PSI( J)-PH| ) ) •$ INF (EN*PS I ( Jl ) 

PHI  -  PHI  ♦  DELTA 

IF  (SPHI  -(PHI  ♦O. 5«0ELT  A  J  I  44*44*43 

44  PHI  »  PHI  -  0. 5»DEL  T A 
ASSIGN  35  TO  JJ 

GO  TO  16 

43  ASSIGN  45  TO  JJ 
GO  TO  18 

45  ZETA  «  ZETA  *2 .0*TAU« S INF ( EN« ( PS  I ( J »-PHl )) «COSF ( EN«PSI ( J M 
ETA  •  ETA  ♦2.0*TAU»SINF(EN*(PSI(J)-PHII)»S1NF(FN»PSI(J)) 

PHI  -  PHI  ♦  DELTA 

ASSIGN  40  TO  JJ 

GO  TO  18 

50  SIPSKJI-  (DELPHI/3. 0)«ALPHA 
S2PSKJ)  -  (  DELPHI  /  3.0 )  »BE  T  A 
S3PSMJ)  -  (  0ELTA/3.0)*ZETA 
S4PSKJI-  (0ELTA/3.0I*ETA 
PSKJ41)  *  PSIIJ)  ♦  DEPSI 
53  CONTINUC 

IF  (KL)  55**0.60 

55  WRITE  OUTPUT  TAPE  2*56 

56  FORMAT (82H0PS I  S1PSI  S2PSI 

l  S3PSI  S4PSII 

00  59  J  -  1 *K*NP 

WRITE  OUTPUT  TAPE  2* 57, PSI ( J I , S IPS  1 1 J » , S2PS I ( J I , S3PSI ( Jl .S4PSI ( J ) 

59  CONTINUE 

57  FORMAT ( 1H  E14.6,4E1 8.6) 

60  I  •  3 

CALL  HARNAL  (PSI. SIPS  I  I 
WRITE  OUTPUT  TAPE  2,61 

61  FORMAT  (  lH0/'//6HQSiPSI//////) 

CALL  PR INCE( L , AQUAY *  BQUAY * AZERO ) 

CALL  HARNAL  (PSI.S2PSI) 

WRITE  OUTPUT  TAPE  2,63 

63  FORMAT! 1H0////6H0S2PSI//////I 

CALL  PRINCE  (L , AQUAY* BQUAY .AZERO I 
CALL  HARNAL  (PSI,S3PSI) 

WRITE  OUTPUT  TAPE  2*65 

65  FORMAT ( IH0////6H0S3PS I//////) 

CALL  PR INCE ( L .AQUAY *  BQUAY .AZERO) 

CALL  HARNAL  (PSI*S4PSI) 

WRITE  OUTPUT  TAPE  2,67 

67  FORMAT ( 1H0////6H054PS I//////I 

CALL  PR INCE ( L , AQUAY,  6GUAY , AZERO I 
80  IM  »  IM  -1 

IF  (IM)  82.62,11 

82  IL  ■  IL  -  l 

IF  (ID  84,84,10 

84  IK  •  IK  -  l 

IF  (IK)  86,86,8 
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86 


I  J  =  I J  -  1 
IT  ||j)  88,88.6 

88  I  I  3  I  I  -  1 

r  (II)  92,92,4 

90  WHITE  OUTPUT  TAPE  2 , 9 1 , PH  I , PS  I ( J  ) 

91  FORMAT  (1H3/////22H  SINGULARITY  AT  PHI  «  E15.6.10H  PSKJI  •  E 15. 61 

GO  TO  eo 

92  CALL  EXIT 

END! i, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0,0, 01 
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SUBROUTINE  PR INCE ( L , AQUAY , BOUAY , AZERO ) 

DIMENSION  AQUAV ( 20 ) • BQUAY ( 20 ) 

WRITE  OUTPUT  TAPE  2,5,AZER0 
5  FORMAT! 1 HO/// 10H0A -ZERO  «  E18.8) 

DO  9  J  =  l  f  L 

WRITE  OUTPUT  TAPE  2 , 7 f J , AQUAY ( J I , J , BQUAY ( J > 

7  FORMAT  I 3H0A- I  2, 3H  -  E18.8.10H  B-I2.3H  =  E18.8) 

9  CONTINUE 

RETURN 

ENO( 1,0,0,0»0*0»0»G,0, 1  ,0,0, 0,0,0) 
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SUBROUTINE  HARNAL  (PSI.FPSI) 

DIMENSION  PSI(36l),FPSI(361)»AQUAY(?0) ,BQUAY(20> 

COMMON  PS  I, F PS  I ,1 ,L,K,DEPS 1, AQUA Y.BQUAY, AZERO 

KN-K-3 

R*1.0 

GO  TO  (5,35,51,1 

5  GRAL  *  FPS  1(1)  ♦  FP  SICK)  +4 .0*  FP S  I  ( K- 1  ) 

DO  12  M*  2  ,  KN  »  2 

12  GRAL  *  GRAL  ♦A.O»FPSI(M)  *2 . 0* FP S  l  ( M* 1 ) 

AZERO  ■  (DEPS I*GRAL )/9. 4247778 
0*0.0 

DO  25  J=1,L 
0*0  ♦l.O 

GRAND*  FPS I ( 1 ) »COSF( Q»PS  I  (  l  )  )  *F PS  I <K I •COSF ( 0* PS  I ( K  )) 

1  ♦4.0*FPSI(K-n#COSF(Q#PSI(K-l)  ) 

DO  20  M* 2 , KN , 2 

20  GRANOGR  ANi)  ♦A .  0  *COSF  ( Q*P  S I  (  M  ))  *F  PS  I  (  M  ) 

1  ♦2.0*FPSI (M*l )«COSF(Q«PS!(M*l I  > 

AQUAY(J)  -  (  (2.0*DEPSI)/9.4247778)«GRAND 
IF  ( 1-3)  25, 22,22 
22  AQU AY ( J ) * A  QUA Y ( JI/2.0 
25  CONTINUE 

IF( 1-3)49,34,34 

34  AZERO  *  AZERO/ 2.0 
R  *0.5 

35  P  *0.0 

00  43  J  * l »  L 
P  *P  41.0 

TGRAL  -  FPSI  m«SINF(P*PSI(l)>  ♦  FPSI  ( K  )  •  S  INF  ( P*P  S 1  (  K  )  ) 
1  ♦4.0»FPSI (K-l ) • S I N F ( P»PSI(K-i  )  ) 

DO  42  M  * 2 • K N, 2 

42  TGRAL  * T GR AL  44.0 »FPS I  (  M )• S INF ( P #P S I ( M ) ) 

1  42.0*FPSI (M4l)*SINF(P»PSI(M4l  )  ) 

BQUAY ( J )  = ( <2.O*0EPSI ) /9 . 424  7778  )*TGRAL »R 

43  CONTINUE 
49  RETURN 

END ( 1 ,0,0, 0,0,0, 0,0,0, 1,0,0, 0,0,0) 
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APPENDIX  D 


Typical  Results  from  Program  of  Appendix  C 
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-.2633196 

01 

.1617436  00 

-.39040)6  00 

.1032606 

ol 

-.1120936 

00 

-.2032326 

01 

.2622006  00 

-.2622006  00 

.196)506 

01 

.0227666 

00 

-.2669906 

01 

.2957606  00 

-.1225116  00 

.2094)96 

01 

.1590316 

01 

-.2200006 

01 

.2719906  00 

-.1469206-06 

.2225296 

01 

.2071416 

01 

-.1506056 

01 

.2103376  00 

.0712426-01 

.2356196 

01 

.2300946 

01 

-.  102799E 

01 

.1)20196  00 

.1320196  00 

.2407096 

01 

.2450526 

01 

-.5570436 

00 

.3043146-01 

.1411146  00 

.2617996 

01 

.2579066 

01 

.7973946-02 

.  U3471E-06 

.1223276  00 

.2740096 

01 

.2522376 

01 

.7633006 

00 

-.3606966-01 

.0901146 -01 

.2079796 

01 

.2132636 

01 

.1551766 

01 

-.5263706-01 

.3263796-01 

.30  1  0696 

01 

.1522306 

01 

.2100716 

01 

-.5194516-01 

.2131646-01 

.3141596 

01 

. 7444096 

00 

.2574796 

01 

-.4131106-01 

•4017O2r-O7 

.3272496 

01 

-  > 1 107066 

00 

.271 095E 

01 

-.2711406-01 

-.1123096-01 

. 3403396 

01 

-.1006026 

01 

.2506006 

01 

-.1416006-01 

-.1416006-01 

.35)4296 

01 

-.1010436 

01 

.219059C 

01 

-.3000306-02 

-. 1209106-01 

.3665196 

01 

-.2453716 

01 

.1620176 

01 

-.1297326-07 

-.0)31246-02 

.3796096 

01 

-.2964116 

01 

.9501416 

00 

.2146046-02 

-.5101006-02 

.3926996 

01 

-.3347366 

01 

.1O1370E 

00 

.3400166-02 

-.3400176-02 

.4057096 

01 

-.3460316 

01 

-.9994796 

00 

.5469646-02 

-.  224502T -02 

.4100796 

01 

-.3140446 

01 

-.2210056 

01 

.0911306-02 

-.1600106-07 

.4319696 

01 

-.2345046 

01 

-.3333306 

01 

.1270096-01 

.5297316-02 

.4450596 

01 

-.1137116 

01 

-.4141996 

01 

.1471536-01 

.1471536-01 

.4501496 

01 

.3530616 

00 

-.4400466 

01 

.1146406-01 

.2767636-01 

.4712)96 

01 

.1941246 

01 

-.4240136 

01 

.90  51  796  -  07 

.4146456-01 

.404)296 

01 

.3395796 

01 

-.3)01006 

01 

-. 2 12603E-01 

.5132726-01 

.4974196 

01 

.4471066 

01 

-.19)0496 

01 

-.5125356-01 

.3125306-01 

.5105096 

01 

.4912226 

01 

-  .4363006- 

-01 

-.0547446-01 

.  3540496-01 

.5235996 

01 

.4501346 

01 

.2030036 

01 

-.1159256  00 

.2007066-06 

.5)66096 

01 

.  3135766 

01 

.  392703E 

01 

-.1320686  00 

- . 547041 E-OI 

.3497796 

01 

.9107916 

00 

.5156026 

01 

-. 1220716  00 

-.1220716  00 

.5620696 

01 

-.1773706 

01 

•329309E 

01 

-.7976676-01 

-.1923726  00 

.5759596 

01 

-.4290576 

01 

.4111016 

01 

-.6943426-06 

-.2460696  00 

.3090496 

01 

-.5073076 

01 

.1773606 

01 

.1104926  00 

-.2667536  00 

.6021306 

01 

-.3090666 

01 

-.1052596 

01 

.2357646  00 

-.2357656  00 

.6152206 

01 

-.4)16746 

01 

-.3420616 

01 

.3510906  00 

-.1454276  00 

.6203106 

01 

-.1762546 

01 

-.4600346 

01 

.  42  7254E  00 

-.1303566-05 

.000 
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S1PSI 


A-ZERO  « 

-•42709232E-01 

A-  1  « 

-.21139412E  00 

B-  1  - 

•2 1096176E  00 

A-  2  « 

-.3711I007E  00 

B-  2  - 

-.U780S91E  01 

A-  3  ■ 

- • 1180 1 338E  01 

B-  3  - 

•27526765E  01 

A-  A  - 

-.21306370E  00 

B-  4  • 

• 19521114E  01 

A-  5  • 

•4902611 5E-01 

B-  5  • 

. 80651021E  00 

A-  6  « 

•87392528E-01 

8-  6  ■ 

•  2631 9832E  00 

A-  7  ■ 

•57828449F-01 

8-  7  • 

• 1 6170567E  00 

A-  8  ■ 

• 17061530E-01 

B-  •  ■ 

•7492B516E-02 

A-  9  - 

•23163780E-01 

6-  9  • 

•97607718E-02 

A-10  - 

. 32263728E-02 

8-10  - 

• 12639331E-01 

S?PSI 


A-ZERO  « 

-• 23103523E  00 

A-  1  * 

-.117315 06 E 

00 

B-  1  ■ 

-.12839404E 

00 

A-  2  = 

•  l 3  1 10791 E 

01 

B-  2  - 

- • 342142  75E 

00 

A-  3  - 

- • 2541 5402E 

01 

CD 

• 

W* 

• 

-. 12022938E 

01 

A-  4  a 

18638599E 

01 

B-  4  • 

-.21291893E 

00 
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A- 

5 

a 

-.773604346  00 

B~ 

5 

9 

•  3 11240906-01 

A- 

6 

a 

-.224941836  00 

B- 

6 

9 

• 772014306-01 

A- 

7 

a 

-• 18027032E  00 

B- 

7 

m 

. 54323212E-01 

A- 

8 

a 

• 56673586E-02 

B- 

8 

9 

. 11827571E-01 

A- 

9 

a 

- .534877 76E-02 

B- 

9 

9 

• 23232483E-01 

10 

a 

-.986976856-02 

B- 

10 

a 

.409252616-02 

S3PSI 


A- 

ZERO  ■ 

-.527499786-02 

A- 

1  - 

-.615971836-02 

8-  1  - 

- .33950984E-01 

A- 

2  « 

.101765426  00 

B-  2  » 

- • 10989538E  00 

A- 

3  - 

.246601706  00 

B-  3  - 

.224819716-04 

A- 

4  « 

•101660756  00 

B-  4  » 

.109860316  00 

A- 

5  « 

-.57524160E-02 

9 

1 

VJI 

n 

.333294816-01 

A- 

6  ■ 

-.527476156-02 

B-  6  ■ 

•  191 36677E-02 

A- 

7  - 

-.406810276-03 

N 

r>- 

i 

ec 

- . 6201 575 3E-03 

A- 

0  - 

.993504626-04 

B-  8  ■ 

-.335058296-04 

A- 

9  ■ 

-.798039946-06 

B-  9  - 

.225071176-04 

A- 

10  - 

-.532949326-05 

B- 10  ■ 

-• 1544 1845E-05 

S4PSI 
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A- 

ZERO  = 

. 191344526-02 

A- 

1  = 

•  32710506E-01 

e-  l 

S 

- .  5346 358 OE -02 

A- 

2  = 

•  1 09828  35C  00 

B-  2 

- 

•  101 56655E  00 

d  - 

3  = 

.224785046-04 

B-  3 

a 

•  246603 1 2E  00 

A- 

4  * 

-.109863426  00 

B-  4 

a 

.101671236  00 

A- 

5  = 

-.333319986-01 

P-  5 

s 

- .  575  3660 5E -02 

A- 

6  = 

- . 19 1 322206-02 

B-  6 

a 

-.527524156-02 

A- 

7  =* 

•  620  123666-03 

R-  7 

2 

-.406560836-03 

A- 

8  * 

.335425  136-04 

3-  8 

2 

.995175566-04 

A- 

9  = 

-.224339046-0'* 

P-  9 

z 

-.620985176-06 

A- 

10  = 

.157688436-05 

B-  10 

3 

-.514121096-05 
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APPENDIX  E 


Program  II  with  Typical  Results 


01  MENS  ION  PS  I (fO),ANG(50),ZETA(6),TW(50) , T (6* 50) ,C0 ( 50 ) *  SI ( 50 ) » 

I A( 20) • 6 1 20) 

COMMON  £L,£MU, AM60A, ETA 

1  READ  INPUT  TAPE  A, 2  *  1 1 , IK, NZ »NM, K, OEPS  I  ,  AO 

2  FORMAT(3I2,2I3«E13.6,E8.3) 

ZN  •  N7. 

CON  ■  6.2831852/ZN 
ZETAID-0.0 
00  10  I  -  2.NZ 
10  ZETAID-ZETAI  I-iUCON 

PSI(l)  -  0.0 
DO  3  J  -  1 .  K 
ANGU)-57.295781*PSt  (J) 

COI  J)-COSF(PSIU) ) 

SIIJ)  -  S INF ( PS  I ( J ) ) 

3  PSI  I  J*1)«PSI  (JUOEPSI 

A  READ  INPUT  TAPE  A , 5 , EMU. AMBOA, I J 

5  FORMAT ( 2E0. 3, 12) 

6  READ  INPUT  TAPE  A, 5, EL, ETA 
BLE*AO» I  EL-ETA) 

WRITE  OUTPUT  TAPE  2 , 1 2 .NZ , EMU. AMBOA , EL ,£TA ,DEPS I 

12  FORMAT  I 1HI  12,1AM  BLAOES  MU  -  F6.3,llM  LAMBDA  -  F6.3.6M  L  - 

1F6.3.8H  ETA  •  F6.3.13H  OELTA  PSI  ■  F8.5/IH0/66H  PSI  TOTA 

2L  WAKE  CONTRIBUTIONS  FROM  INDIVIDUAL  BLAOES/1HO) 

DO  30  J  •  l.K 
T  W ( J )  •  0.0 
00  25  I  -  l.NZ 
PZA-PS1  (  JUZETA(I) 

I F ( I  —  1 )  13,13,16 

13  TMNU-FMU»CO( Jl 
THDE«EL*EMU*SI(J) 

IF(ABSFITHOE)-. 00001 )  20,1A,1A 

1A  OELTA-ABSFIATANHTHNU/THDE)  ) 

IF(THDE)  15,115,115 

15  0ELTA-3.1A15926-DELTA 

115  IF(THNU)  116,117,117 

116  OELTA  -  -  OELTA 

117  T( I , J  1  -  * 1.0-SINF (OELTA) ) /( ( EL-E T A) •COSF  I  DE LT A )  ) 

GO  TO  22 

16  CALL  0UD1(SI(J),PSI< J),PZA,KER,PHI) 

PZP-PZA-PHl 

I F ( KER-2  )  17,60,17 

60  WRITE  OUTPUT  TAPE  2,61,ZETAII) 

GO  TO  20 

61  FORMAT ( 25HONO  CONVFRGFNCE  AT  Z  FT  A  -  F 1 0. 7  » 

17  D  -  EMU*PZP 

Z  •  |AMB0A*PZP-BLE)««2 
ERN  •  EL*»2-(0*SIU)  )**2 
IF(ERN)  20,21,21 

20  TII.JI  •  0.0 
GO  TO  22 

21  ERN  ■  SORTFIERN) 

YC«0*C0(JUERN-ETA 
THNU*-EMU»CO$F (PHI ) 

THDE-FL*EMU»SINF(PHl ) 

IF(ABSF(THOE)-. 00001 )  20,70,70 
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70  PS0«ABSF(ATANF( THNU/THDE  I ) 

I F ( THOE )  71,72,72 

71  PS0»3.1415926-PSD 

72  IF(THNU)  73,74,74 

73  PSO--PSD 

74  0ELTA«PHI-P$!  (  JJ-P50 
UTT«YE«COSF!OELTA| 

T(  I,J)*2.0*UTT/(Z^UTT»*2) 

22  TW(JI*T(I, J)4TWC J) 

25  CONTINUE 

30  WRITE  OUTPUT  TAPE  2,33, ANG(  J ),TW(J),(T(I,J),I»l,NZ) 

33  FORMAT ( 1H  F8*1»F13*5»F15*5»5F12*5) 

CALL  HANEWIPSI  ,TW,K,NH,0EPSI ,AZERO, A,B) 

WRITE  OUTPUT  TAPE  2, 44 , AZERO , ( N, A( N ) ,N ,B(N ) , N- l , NH I 
44  FORMAT! 18H1HARM0NIC  ANALYS IS/1H0/9H0A-ZER0  -E 1 2. 3/ I 3M0A- 1  2 , 1M*  E 
112*5, 5X,3H  B- 1 2, IH«  E12.5U 

PUNCH  101, AZERO,  (A(J),J«1,  NH) ,  (BUI,  I-1,NH) 

PUNCH  100, CTW! J),J-1,K) 

100  FORMAT (6E1 2 • 5 ) 

101  FORMAT ( 5E14*  7  I 
IJ  ■  IJ  -  1 

I F ( I J )  53,53,6 
51  IK  •  IK  -  1 

IF! IK)  54,54,4 
54  II  «  II  -  l 

IF(II)  56,56,1 
56  CALL  EXIT 

END! 1 ,0,0, 0,0, 0,0, 0,0, 1,0,0,0,0,01 
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SUBROUTINE  OUDl < S I , PS  1 , PZ A,KER , PH  I ) 

COMMON  EL*FMU 

NC  »  l 

KFR  *  l 

Cl-EMU»$l 

PHI  -  PSI 

2  FFOD-C1 •< PHI -PZ A  )*EL*SINF{ PH I-PS I ) 

EFPOO*CUFt*COSF  (PH  I -PS  I) 

AO  -  EFOD/EFPOD 
I  F ( ABSF ( AO)-. OOOOI )  10, 10, A 
4  NC  ■  NC  ♦  1 

I F ( NC“ 10 1  6,8,6 
6  KFR  -  2 

GO  TO  10 

8  PHI  «  PHI  -  AO 

GO  TO  2 
10  RETURN 

E  NO ( 1 ,0,0, 0,0,0, 0,0,0, 1 ,0,0,0,0,01 
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SUBROUTINE  HANEWIPSI ,T,K,NH,OEPSI, AZERO, A,B) 
01  MENS  ION  PSI(50),T(50)«A(20),B(20),W(50) 

0  ■  OE PS  1/9*4247770 
AZERO  «  T(1)+T(K)+4«0*T(2) 

00  2  II  -  1,NH 
A(II)  -0*0 
2  S(II)  >  0.0 

Min-i.o 

W(K)-1.0 

W<2)«4.0 

KKK  -  K  -  l 

DO  4  JJ  •  3  f  KKK  , 2 

W( JJJ-2.0 

W( JJ*1 )-4*0 

AZERO  •  AZERO  ♦2.0*T ( JJ ) ♦4.0*H JJ*l) 

4  CONTINUE 

AZERO  •  0.5*Q*AZER0 
00  8  II  -  1 1 NH 
ELL  -  II 
00  10  NB  «  l,K 

A( II )-A( lI)*W(NB)*TCNB)*CO$FCFU*P$I(NB) ) 

10  B(II)*B{II)*W(NB!"7(N8)«SINF{FIL*PS1(N&)) 

A ( 1 1 ) «Q*A ( II) 

B(II>  *  Q*B( 1 1 1 
8  CONTINUE 
RETURN 

ENO( 1,0, 0,0, 0.0, 0» 0,0, 1,0, 0» Ot 0,0) 
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4  SHOES  MU  •  .200  LAMBDA  • 


029  L  ■  1.000  ETA  -  .950  DELTA  PSI  •  .26180 


PSI 

TOTAL  WARE 

CONTRIBUTIONS  FROM 

1N0IVI0UAL  BLADES 

0. 

26.93724 

16.39608 

3.54507 

2.. 96902 

2.02700 

15.0 

20.92674 

16.66103 

6.08904 

3.37927 

2.39660 

30.0 

31.75953 

17.09723 

7.15731 

4.24545 

3.25956 

45.0 

37.94855 

17.67494 

9.07674 

5.96033 

5.22854 

60.0 

51.02715 

18. 361 79 

12.79333 

9.09772 

9.96032 

75. 0 

45.90260 

19.15115 

21.52012 

15.14319 

-9.91170 

90.0 

24.75162 

20.00000 

27.70366 

-15.42452 

-7.52752 

105.0 

-22.80051 

20.08647 

-29.71257 

-0.97072 

-4.99569 

120.0 

-3.02990 

21.77725 

-16.90221 

•6.17026 

-3.73468 

135.0 

2.99334 

22.63091 

-11.01830 

-4.00010 

-3.01116 

150.0 

7.34065 

23.39960 

-9.41402 

-4.06156 

-2.57856 

163.0 

9.05217 

24.00812 

-0.16925 

-3.64970 

-2.33700 

100. 0 

11.10604 

24.39608 

-7.57625 

-3.46602 

-2.24696 

195.0 

11.22926 

24.48546 

-7.46246 

-3.47990 

-2.31305 

210.0 

10.03049 

24.21600 

-7.02017 

-3.71366 

-2.62360 

225.0 

6.01905 

23.56301 

-8.00556 

-4.25384 

-3.60455 

CONVERGENCE  AT  ZETA  * 

4.7123889 

240.0 

6.36070 

22.96473 

-10.91124 

-5.28400 

0. 

255.0 

16.09022 

21.32429 

-15.67933 

-7.02600 

10.27926 

CONVERGENCE  AT  ZETA  ■ 

4.7123089 

270. 0 

-0.32343 

20.00000 

-19.19659 

-9.12088 

0. 

285.0 

14.77329 

10.75796 

20.03869 

-9.07433 

-14.94902 

CONVERGENCE  AT  ZETA  • 

4.7123889 

300.0 

30.60802 

17.72678 

10.64250 

10.23955 

0. 

315.0 

35.21124 

16.97518 

7.17173 

5.05356 

6.01077 

330.0 

20.30611 

16.51800 

5.84963 

3.37712 

2.64136 

345.0 

26.71420 

16.33622 

5.43521 

2.92165 

2.02112 

360.0 

26.93724 

16.39608 

5.34507 

2.96902 

2.02708 
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HARMONIC  analysis 


A- 

ZERO 

*  • 1 8087E  02 

A- 

l* 

.159470 

02 

B-  1« 

•91212E 

00 

A- 

2* 

.29015E 

01 

8-  2- 

.62804E 

01 

A  - 

3*  - 

•  1 1 5500 

02 

8-  3- 

-.11793F 

01 

A- 

A*  - 

•  3591 1 E 

01 

8-  4* 

- .73613E 

01 

A  _ 

5* 

•44264E 

01 

B-  5* 

• 16926E 

01 

•\- 

6« 

• 19471E 

01 

8-  6* 

• 74949E 

01 

A- 

7-  - 

• 18754E 

01 

8-  7- 

-.19209E 

01 

A- 

8* 

•  28 1 1  IE 

00 

R-  8- 

-.56557E 

01 

A- 

9« 

.349810 

01 

B-  9* 

•24231E 

01 

A- 

10*  - 

•  1 2456F 

01 

8-10* 

•34025E 

01 
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NO 

CONVERGENCE  AT  ZETA  - 

3.1415926 

300.0  7.77835 

-2.95908 

0. 

0. 

NO 

CONVERGENCE  AT  2ETA  - 

3.1415926 

NU 

CONVERGENCE  AT  ZETA  ■ 

4.7123889 

315.0  -13.50982 

-2.72228 

-10.78755 

0. 

NO 

CONVERGENCE  AT  ZETA  ■ 

4.7123889 

330.0  20.00500 

-2.62685 

14.51674 

6.11511 

345.0  20.60732 

-2.63173 

15.87829 

4.73445 

360.0  23.23921 

-2.70813 

17.66847 

5.29967 

10.73743 


0. 


0. 

2.62631 

2.97920 
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HARMONIC  ANALYSIS 


A-ZERO  -  - • 73798E  01 


A-  1- 

.97773F 

01 

B-  1- 

-  .9134  2E 

o: 

A-  2« 

•  7391 5E 

01 

B-  2- 

.  <f&?76E 

01 

A-  3- 

•  1001 2E 

02 

B-  3« 

. 40O65E 

00 

A- 

• 76078E 

01 

B-  4* 

.3606  3!'. 

on 

A-  5- 

• 30506F 

01 

B-  5- 

•  167 I 2E 

01 

A-  6- 

•47120E 

CO 

B-  6* 

• 13112E 

01 

A-  7* 

38262E 

01 

B-  7- 

•  1 3 1 26E 

01 

A-  P- 

-• 63643E 

01 

B-  8* 

•59635F 

00 

A  —  9* 

-• 53324E 

01 

B-  9* 

- • 40838E 

00 

A-10- 

-• 30983E 

01 

B-10- 

-.61389E 

00 
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APPENDIX  F 


Program  III  with  Typical  Results 


OIMENSION  Et  MIOI  ,  THETA  |  10)  ,FI10),GU0)  ,PS1  ISO)  , AND  I  50)  ,SI  (50)  ,C0 
l(50),Cr(50) tST(50) , C  T  T I  50 ) ,STT(50),C(21 )  >0(50) «  A  < 1 j  ,  1 0 . 2 1 ) , B ( 1  J ,  10 
2 ,50)  ,0(20,?  I  )  ,  Z  (  20)  ,  T  IMPS  (  10)  , EL AMI  10,  2  I  )  ,PL  (  10, 50)  ,GC  f  1 6  )  ,GS(  10  )  , 
3FLC ( 8  )  , ELS ( H  ) 

1  READ  INPUT  TAPE  A  ,  2  ,  DE  PS  I  ,  BOER  ,  t  ANMUi  EMUi  C?  I  1  .  Mi  KP  ,NH  ,NGP  I  ?P 

2  F0RMAT1E1 3.8.4E8. 3,613) 

N2*  NH/2 

MN  *  M  ♦ 1 
ROB  *  1.0/80ER 
DEMU  *  0 . 5»  EMU 
E  2M  *0 . 5  * EMU* EMU 

3  READ  INPUT  TAPE  A , 4 , ( E T A (  I  )  ,  THE T A ( 1  ) , F (  I  > , G ( I  ) , I  =  1 , M ) 

4  FORMAT  ( Af  8 . 3 ) 

PSI  (  1  )-0.U 

ID  *  d 

DO  6  J  1  1  *  KP 

ANG(  j)*57.20,:’Bl«P$I(j) 

C0( J )  »  COSFIPSI ( J)  1 
SI ( J)  «  SINFIPSII J)  ) 

CTI J)  «  COSF ( 2 . 0»PS I ( J ) ) 

ST  IJ)  *  SINF(2.0«PSU  J)  ) 

CTTI JI*C0SF(3.0*PSI(J)  ) 

PSI(JM)  =  PSIIJ)  ♦  OEPSI 

6  ST  T I J  )  =  S 1 NF I  3 . 0 *P S 1 1 J  )  ) 

DO  7  I J I  *  1  a  M 

7  TIMES!  I J  I  )  =  F T A  I  I J I  ) • THE T A  I  I  J 1 ) 

00  20  I  *  1,«  ‘ 

DO  13  J  >  1 » MN 

READ  INPUT  TAP';  A  ,  8 , I  A (  I , J , K  ) , K *1 , NH ) 

READ  INPUT  TAPc  A , 9 , ( B l I , J , K  I , K *  1 , KP ) 

8  FORMAT ( 5E 1  A . 7 ) 

9  FORMAT  (6E12.5) 

GO  TO  ( 13, 10) , I2P 

10  READ  INPUT  TAPE  A,8, (C(l) »  L  3  1 ,  NH  ) 

READ  INPUT  TAPE  A  ,  9  ,  I  D  l  t. )  ,  L  *  1 ,  KP  ) 

DO  11  K  *  1 , NH 

11  A ( I . J . K I  *  At  I,  J,K)  ♦  C  (  K  I 
00  12  L  *  1 , KP 

12  9(T,J7U~=  MI.J.U  ♦  OUT  ~ 

13  CONTINUE 

Q< I,MN)*ROD*f ANMU 
DO  18  N  =  1 , M 
00  15  K  *  1 , NH 

15  AU.N.K)  =  AII.N.K)  -  \I  I  ,  l  ,K  ) 

DO  16  i  *  1 ,  K  P 

16  BII.N.L)*  BCI,N,L)-B(l,N*l,L> 

01  I iNl'AI I ,N, 1  I 

1FII-N)  18,17,18 
17  01  I , N) -  0 ( I , N )  ♦  ROB 

18  01  I ,MN)*0< I ,MN ) *AC 1 ,N, l )• T IMES(N) 

20  CONTINUE 

CALL  CROUT I J J , M , 0 , Z I 
1FIJJI  2A, 24,22 
2?  WRI TE  OUTPUT  TAPE  2,23 

23  FORMAT ( 38H1 NO  SOLUT I ON--D I  AGONAL  ELEMENT  -  7FR0.) 

GO  TO  64 
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24  WRITE  OUTPUT  TAPE  2 , 25 , I E T A (  I) , I  * l , M ) 

25  FORMAT! 28HILOAD  COEFFICIENTS  FOR  F  T  A  = , 9 ( T 6 . 3 , 1 H , ) / / ) 

WRITE  OUTPUT  TAPE  2,26,lD,(Z(l)*I-l,M) 

26  FURMAT! I3HOHAKMONIC  NO.  ,12,9^13.5) 

DO  28  K  *  1  ,  M 

28  TIMES!K)=HOCR*t TIMES! K)-Z!K)  ) 

DO  35  J  =  2 *  NH 
ID  =  J  -  1 
DO  30  I  =  1  •  M 
ELAM! I t J)  =  O.U 
DO  30  K  =  1  *  M 

30  ELAM(  I , J )  =  E l A  M l I,J)+T|MES!K)«A!  I  ,  K  ,  J  ) 

35  WRITE  OUTPUT  TAPE  2 , 26  ,  I  U ,  !  EL  A  V  \  N  ,  J  )  ,  N  =  l~,  M  ) 

DO  41  I  =  1  ,M 

WRITE  OUTPUT  TAPE  ?,42,ETA(I) 

42  FORMAT  I  22H1  LAMBDA!  PS  I)  FOR  ETA  =■  F6.3//I 

DU  41  J  =  1 ,KP 

PLII.J)  =  T  AN MU 
DU  43  K  =  1 , M 

43  PL!I,J)=PL!l,J)*TlMrS!K)*RlI,K,J) 

WRITE  OUTPUT  TAPE  2 , 44 , ANG ! J ) , PL ( I , JT 

44  FORMAT ! 1H  F7.?,l l 3 . 6 ) 

41  CJMTI'JUF 

WRITE  (JU  T  PU  T  TAPE  2,45 

4  5  FORMAT  I  4H1FTA,9X,  4HL3-C,  1 1 X , 4HL 3- S , i  l  X , 4HL4-C  ,  1 1.  X  ,  4HL4- S  ,  1 1  X  ,4HL5 

1-C , 1 1 X, 4HL5-S, 1 1 X, 4HL6-C, l 1X.4HL6-S/  ) 

DO  50  J  =  l,M 
DO  46  I  =  2,7 
NNN  =  N?  ♦  I 

GC! I )=F! J)*ELAMIJ, I ♦ 1 ) -G ! J ) »EL AM ! J , NNN ) 

46  GS  (  I  )=FIJ)*ELAM!J,iMNN)  *G  !  J  )  *  EL  AM  {  J  ]  I  ♦  1  j 
DO  47  K  =  3,6 

FLClK)  =  CF*l  FT  A  (  J  1  «GC  (  K  )  ♦DEMU*  (  GS  ( K  ♦  lT-G$  ( K-  1  )  )  ) 

47  ELSIK )  =  C  F  •  ( FT  A! J) • G  S I K ) ♦OEMU* l GC t K- 1 ) -GC IK* 1 ) )  ) 

WRITE  OUTPUT  fAPE  2 , 49  ,  F T A  {  J )  ,  ( ELC  !  L  )  ,  EL  $ « L ) , L  =  3 , 6 ) 

49  FORMAT !  1H0F6. 3, HE  1 5.6) 

50  CONTINUE 

52  READ  INPUT  TAPE  4 , 54 , A Z FRO , AON E , BONE 
54  FORMAT! 3E8. 3) 

DO  63  I  =  1 ,M 

WRITE  OUTPUT  TAPE  2, 56, ETA!  I  j 

56  FORMAT  !  6H 1  C  T  A  =F6.3/4H0PSI  ,  l  IX,  <*HL  I  I  )  ,  12X,4HL  (  G  )  ,  1  OX  ,  9HL  tl)  ♦  L  I  G  )  , 

15X,26HL(GM  ) -HARMON  ICS  EXTRAC  TF  0/  /  ) 

C  =  E  r A!  I ) • *2 

Cl* THETA! I >*<C*E2M)-ETA{ I )*FMIJ*A0NE 
C 2*2.0* EMU* THETA! I ) *ETA! I )-AONE*!C» l.5*F2M) 

C3  =  0(JNlE»!C^.5»E2M)-EMU*AZERO»f:l  A!  I  ) 

C4  =  EMU«R0Nt *ET A f  !  )-AZER0*E2M 
C 5*  E MU • AONE • E  T  A !  I  )-THE T A l  I  )  *L 
C  6=  0 . 5  *E  ?M • AUNT 
C  7s  -0.5*E2M*PUNE 
B1*ETA(I)*Z!I MDTMU-ELAM! I,N?*1) 

P?  =  F  T  A !  I  ) • T  L  AM (  1»Z)*DEMU*FLAM!  I , N2 ♦ 2 ) 

B3=FTA{ I >*ELAM( I , N2 ♦ l ) ♦ EMU* Z ( I ) -OE MU* C L AM { 1,3) 

B  4  =  F  T  A!  I  )*CLAM(  1,  3)-DEvU*(ELAM(  I ,  'J  2  ♦  I ) - E  L  A  M ( I  ,  N  2  ♦  3  ) ) 

B5*ETA{ I )*ELAM( I,N2*2)*DEMU*(ELAM! I,2)-ELAM( 1 ,4) ) 
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00  63  J  =  l  ,  K  P 

'l  =  CF  •  ( E  r  A(  I)  ♦EMU»SI  (  J  )  )  _  ___  _ 

ELCYE=-K«PL ( 1 , J ) 

FLG=CF*(Cl*C2*SI(J)*C3*CO(J)  *C4*  S  T  {  J  1  *C5»CT  (  J  1  ♦C6«STT  (JJ*C7*C1  T  l  J) 
1  ) 

SUML  =  FLG  +  CL  EY 6 

HEL=  +  ELEYF  +  CF»(Bl-»-B2»CO(J)4-93*SI(J)  +B4«CT  IJ1+B5«ST(J)  1 
WRITE  OUTPUT  TAPC  2 , 58 , ANG ( J ) , EL F YF , ELG, SIWL ,HF L 
58  FORMAT ( 1 H  F7.2,4E16.6) 

63  CONTINUE 
66  NGP  =  NGP  -  I 

IF(NGP)  66,66,52 
66  11=11-1 

I F {  I  I  )  68,68, 1 
68  CALL  EXIT 

ENO ( 1,3, 0,3, 0,0, 0,0,0, 1,0,0,0,0,01 
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SUBROUT  I  NC  CROUM  JJ,M,  A,  Z  > 
DIMENSION  A ( 20 *  2 1 ),Z(20> 

J  J*  -2 
MN  -  M  M 
DO  27  I  *  1 » M 
NN  =  l 

00  25  J  =  l,MN 
SUMA  =  0.0 
IF(J-l)  25,2*5,9 
9  IF  ( 1-1 )  19,19, 12 

12  DO  15  K  =  l , NN 

SUMA  =  SUMA-f  A(  I  ,K)  *A(K,J) 

15  CONTINUE 

17  A(  I  ,J)=A(  I,  Jl-SUMA 

IF  (l-JJ  19,26,26 
19  IF  (All, I))  23,21,23 

21  J  J  =  A 

CO  TO  35 

23  A ( I, J)  =  A (  I , J  I / A ( 1,1) 

26  I F  C  I -1 -NN I  25,25,  24 

24  NN  =  NN  *1 

25  CONTINUE 

27  CONTINUE 
IN  =  M 

ZI  IN)  =  A(M,MN) 

28  IN  =  IN  -1 
IF(IN-l)  35,29,29 

29  SUMB  =  0.0 
KK  =  IN  ♦  1 

DO  31  K  =  KK , M 

SUMO  =  SUMB  ♦  Z<  K ) * A (  IN,K  ) 

31  CONTIMUE 

Z ( IN)  =  A ( IN , MN )  -  SUMB 
GO  TO  28 
35  RETURN 

END ( 1 ,0,0,0, 0,0, 0,0,0, 1 ,0,0, 0,0,0) 
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LOAD  COFFF I C I  ENT  S  FOR  ETA  = 

.950,  .850, 

.  7 oO,  .600 

t  .300, 

HARMONIC 

NO.  0 

.  3274  7E-01 

.26967E-01 

.242870-01 

. 26822E-0! 

.20601 6-01 

HARMONIC 

NO.  I 

.222310- )1 

.218956-01 

. 16365F-01 

-.  lP054E-dTT 

-.21 749E-02 

HARMONIC 

NO.  2 

.414600-32 

-. 44885E-02 

-.144190-01 

-.43768F-07 

- . 987S5r-03 

HARMONIC 

NO.  3 

-.132050-01 

-. 92275 E -02 

- • 38956E-02 

-. 42504E-02 

-. 45973E-02 

HARMONIC 

NO.  4 

-.88347  f-02 

- .  24R48E-0  3 

- .24182E-02 

. l5*70F-0? 

-.712726-03 

HARMONIC 

NO.  5 

.254960-02 

- . 65 1 05E -  03 

-.268470-02 

.21 982E-02 

-.29281 E-0? 

HARMONIC 

NO.  6 

•25599E-J2 

.41512E-03 

• 26852c-02~ 

~'Z 100516-52' 

-.3281 7E-0? 

HARMONIC 

no.  r 

-.2291 5E-J? 

. 12095E-02 

- .945350-03 

-. 11379E-02 

— . 9  2 1 76F-03 

HARMONIC 

NO.  9 

-.22 4690-03 

.48575E-05 

.144196-02“ 

“-:39038E-02“ 

“ - . 204 79E-0? 

HARMONIC 

NO.  9 

.332020-02 

- . 696  7  9E - 04 

-.151426-02“ 

-.436&4E-03 

"-;i4484E-0? 

HARMONIC 

NO.  10 

-  .37242E-02 

. 14563E-03 

.389016-05 

-.1 7468E-C2 

-.T9227E^02“ 

HARMONIC 

NO.  1  1 

.  16680  E- 02 

-.  78  30  2E -0  3 

-.668166-03' 

““.I05i?C-02“ 

‘-:*9356E-03  “ 

HARMONIC 

NO. 12 

-. 12904E- 03 

-.22413C-02 

.20340E-02 

29061 6-02 

'“.263156-03 

HARMONIC 

NO. 13 

- . 94060  E-0  3 

.21760E-02 

-,24570E^4T 

-.  13821E-03 

-.131916-0? 

HARMONIC 

NO. 14 

•24553E-J3 

-.184306-02 

-  . 10530E-02 

.544966-04“ 

. 1 l 778E-02 

HARMONIC 

NO.  15 

•91562E-U4 

-.208826-02 

.  40499E-02 

-.681986-03 

-. 12945E-0? 

HARMONIC 

NO.  16 

-.1 3683E-32 

.324386-02 

-  . 3004  IE-02 

. 14621E-02 

- . 49949E-0  3 

HARMONIC 

NO.  17 

-. 10085E-0? 

. 31465E-04 

-.  13381E-0? 

•41457c-03 

. 86297F-03  ( 

HARMONIC 

NO.  18 

.464  77  E- 03 

-.239846-02 

.22857E-02 

-.  16059E-0? 

-.  1  0907  E-0? 

HARMONIC 

NO.  19 

-.  1574  1  E-(5  3 

.221886-^2“ 

-  .  24689E-02 

-.601 96E-05 

.627706-35 

75 


L  4P3D  A  ( PSI )  F HR  F J A 


950 


0. 

7  .  d0 
15.00 

22. 50 
30.00 
37. 5C 
45.00 

52.50 
60.00 

67.50 
75.00 

82.50 
90.00 

97.50 
105.00 

112.50 
120.00 
127. dO 
135.00 

142.50 
150.00 

157.50 
165.00 

172.50 
180.00 

187.50 
195.00 

202.50 
210.00 

217.50 
225.00 

232.50 
240.00 

247.50 
255.00 
262. oO 
270.00 

2  77. jO 
285.00 

292. 50 
500.00 

307.50 
315.00 

322.50 
330.00 

3  3  7.30 
345.00 

352.50 
360.00 


.  3  9  300  5E  -  0  l 
.4 12342E-01 
.401 294E-C 1 
•  4613t<7E-Ol 
. 575C40C-01 
•629770F-01 
. 627435 E -01 
.62?  757E-0 1 
.5H1023E-C1 
. 51 i076f-01 
. 397935E-CI 
. 324955E-U1 
. 350328E-01 
-. 1 05605C-01 
-. 1 21480E-0? 
.653  595E-02 
.992  324E-0? 
. 1 19370E-01 
. 135129F-01 
. 148359E-U1 
. 160976E-01 
.  1  7  !  5  8  4  E  -  G  1 
. 180841E-01 
. 188912E-01 
. 1 95956E-0 I 
. 202024F-0 1 
.207083E-01 
.21 1029C-01 
. 2 1 366 1 E -0 1 
.21468  IE-01 
.213633E-01 
.209663E-01 
. 200970E-G] 
. 18391 7E-01 
. 151 204E-01 
. 629840E-O2 
. 1 19347E-0? 
.439814E-01 
.456505E-01 
. 549819E-01 
.609663E-01 
.65887 3E -01 
.6  77144E-01 
.61855 3E-01 
. 540788E-01 
.  4  3962  7E  1 
.4  76748E-01 
.454590E-01 
. 393006E-01 
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LAMBDA ( PS  I )  FOR  FTA 


8  30 


0. 

7.50 

15.00 

22.50 
30.00 

37.50 
45.00 

52.50 
60.00 

67.50 
75.00 

82.50 
90.00 

97.50 
105.00 
L12.S0 
120.00 

127.50 
135.00 

142.50 
150.00 

157.50 
165.00 

172.50 
180.00 

167.50 
195.00 

202.50 
210.00 

217.50 
225.00 
232.  *>0 
240.00 
247.^0 
255.00 

262.50 
270.00 

277.50 
286.00 

292.50 
300.00 

307.50 
315.00 

322.50 
330.00 

337.50 
345.00 
352.30 
360.00 


. 32  7020E-01 
. 36698  3E-0 1 
•377216E-01 

•  3  54  52  4  E -0 1 
. 37427 3E -01 
.4  76039E-0 1 
.484284E-01 
. 5  07238E-01 
. 5 1 9439E -0 l 
.465539E-01 
. 4  09969E-01 
. 345839E-01 
. 280900E-01 
. 194914E-01 
.2  J109JE-01 
. 340384E-01 
•693205E-0? 
.28477  IE-02 
.923531E-03 
. 14581 7E-03 
. 20  745  3E -02 
.4  96284E-02 
.72  7524E-02 
. 8  7453  IE-02 
.960222E-02 
. 10054  3E -0 1 
.101577E-01 
.99172  7E-02 

•  94599  3F -02 
.903670E-02 
. 86239 0E-02 
. 754002E-02 
. 5  34030E -02 
.2O0572E-02 

-.264776E-02 
.  33588  3E--01 
. 36482  3 E  —  0 1 
.34231SE-01 
.392416E-01 
.489399C-01 
.504342E-01 
.496072E-01 
. 50320HE-01 
. 512479E-01 
.45652  8E-01 
.405146E-01 
. 399595E-01 
.37167 3 E-Ol 
.327020E-01 
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LAMBDA ( PS  I  )  FUR  ETA 


750 


0. 

7.50 

15.00 

22.50 
30.00 

37.50 
4  5.00 

52.50 
60.  u0 
67.60 
75.00 

82.50 
90.00 

97.50 
105.00 

112.50 
120.00 

127.50 
135.00 

142.50 
150.00 

157.50 
165.00 

172.50 
180.00 

187.50 
195.00 

202.50 
210.00 

217.50 
225.00 

232.50 
240.00 

247.50 
255.00 

262.50 
270.00 

277.50 
285.00 

292.50 
300.00 

307.50 
315.00 

322.50 
330.00 

337.50 
345.00 

352.50 
360.00 


. 1H2S1 5E-0* 
.  23  1  35  3E- .1 
.21903  5E-01 
.31131  >E-01 
. 3  55044F-0 1 
. 37?45oE-01 
•464060L-01 
.44017 /E -0 1 
.4 l 3394 F -01 
.419334E-01 
.4  33239E-01 
.3  1174  iE-Ol 
. 334640E-01 
. 3l748oE-0  I 
.21099 OE-Ol 
. 209385E-0 l 
.2J7185E-01 
•  240902E-0  1 
.23619 OE-Ol 
.  19972  IE-01 

-.940046E-02 

-. 155745E-01 
. 966265E-03 
•250352C-0? 
.351 395E-03 
. 220  f 59E-02 
. 1 75090E-0? 

-.391 148E-0? 

- . 3  96602  E -02 
.305269E-0? 
. 146388E-02 
.760257E-0? 
.1806610-01 
.421 307E-01 
. 36971 5E-01 
. 333526E-01 
.  32  1  32  7E-0 1 
.33H91 IE-01 
. 4  44184E  —  01 
.4545L’OE-01 
.42221 7E-01 
.467704E-01 
.498476E-01 
. 3621 36E-01 
.263378C-01 
. 24509  3E -0 1 
.239612E-01 
.21 3221E-01 
.  1  82  5 1 5  E -0 1 
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LAMBDA ( PS  I )  F OR  ETA 


6  J  J 


0. 

7.50 

15.00 

22.50 
30.00 

37.50 
*5.00 

52.50 
60.00 

67.50 
75.00 

82.50 
90.00 

97.50 
105.00 

112.50 
120.00 

127.50 
135.00 

142.50 
150.00 

157.50 
165.00 

172.50 
180.00 

187.50 
195.00 

202.50 
210.00 

217.50 
225.00 

232.50 
240.00 

247.50 
255.00 

262.50 
270.00 

277.50 
285.00 

292.50 
300.00 

307.50 
315.00 

322.50 
330.00 

337.50 
345.00 

352.50 
360.00 


. 145676E-01 
. 1 84266F-01 
. 2  3  3498F-0 1 
•230004F-01 
•207H68F-01 
.  18764  3E-01 
. 155022E-01 
.2325  1  IE-01 
•478714F-01 
.423  764E-01 
. 34 2  30  3E -0 l 
. 334754E-01 
. 320067C-0 1 
.24161  IE-01 
. 2  76462E-0 1 
.304692E-01 
. 2  734H8C-0 l 
. 164054E-01 
. 1 1 7842E-0 1 
. 1 77455E-ol 
.220762E-01 
.26171  fZ-Ol 
.2 47810E-01 
.246358E-01 
•251G88E-01 

•  266  3?  7E-0 I 
. 3L0334E-01 
. 3 l 1 l 55E-0 l 
.286345E-01 
•277257E-0! 

•  ?804i)  1E-0  l 
.290208E-CI 
.2P5897E-C I 
.256500E-01 
.251 185E-01 
.259390E-01 
. 306452E-0 l 
.349176E-01 
•362258E-01 
. 3  71660E-01 
.386l04r-0! 
.299190E-01 
.230312E-01 
.24  7650E-U 1 
.262 140E-01 
.2  70 83  / E  —  0  1 
.  28 1067E-Gl 
.2113?  7 E  —  0 1 
.  1456  7oE-0  1 
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LAMBDA ( PS  I)  FOR  ETA 


3  jO 


0. 

7.50 

15.00 

22.50 
30.00 

37.50 
45.00 

52. 50 
60.00 

67.50 
75.00 

82.50 
90.00 

97.50 
105. uO 

112.50 
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